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1. INTRODUCTION 


The principle emphasis in this project has been upon 
the following specific tasks: 

(1) The design and calibration of a neutron 
monitor to measure the cosmic-rav neutron 
leakage from the earth's atmosphere at 
altitudes from 400 to_1000 km over a range 
of latitudes from 0<X<90°. 

(2) To determine latitude dependence of the 
cosmic-ray leakage flux. 

(3) Measure the altitude dependence of the 
cosmic ray neutron leakage flux. 

(4) Measure the energy spectrum of the cosmic- 
ray neutron leakage flux in the energy range 
from .01 to 10 MeV. 

(5) To search for solar neutrons during both 
quiet times and solar flare events. 

(6) To evaluate the neutron flux from the inter- 
action of solar protons in the polar 
latitudes. 

(7) To determine the angular distribution of the 
neutron leakage flux as deduced by measure- 
ments of the altitude dependence. 

(8) Verify the solar modulation of the cosmic-ray 
source for the neutron leakage. 


There has been no attempt yet to study the folloxving 
problems : 

(1) The magnitude and fluctuation of energetic 
charged particles over the poles. 

Since the detector on OGO-VI had charged- 
particle detectors with a very large geometrical 
factor, we can determine the fluctuations in the 
low-energy charged-particle fluxes over both 
poles. Such fluctuations are associated with 
geomagnetic cutoff perturbations and solar proton 
fluxes. So far no attempt has been made to 
unravel these affects. 
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(2) The anisotropies in the protons above about 
70 MeV. This can be measured using the 
coincidence banks in the charged particle 
detector array. However, no studies have 
been made. Again the results should be 
rather interesting because our detector 
efficiency is large. We have observed 
charged-particle anisotropies when the 
detector passes into the South Atlantic 
magnetic anomaly. Only a preliminary in- 
vestigation has been made. This was done 
because these data from the anomaly must be 
eliminated from our analyses of the cosmic 
ray neutron leakage flux measurements. During 
such passes there is a large amount of neutron 
production affects by the energetic protons trap- 
ped at these altitudes. 

(3) Short duration counting-rate bursts. These 
might be present in the charged-particle data 
from solar x-ray fluxes incident on the 
detector. This is a large task and involves 
reworking the primary data tapes. It does not, 
however, seem warranted at this point. 

Since most of the proposed analyses for this project 
have been published, we have included several reprints in 
this report. Other results are already in preprint form, 
copies of which are attached, or are contained in a Ph.D. 
dissertation by S. Ifedili. We have not included any 
discussion of the details of the design criteria, the methods 
of construction of the neutron detector, nor of the 
calibration at various particle-accelerator facilities and 
with radioactive sources. These details can be found in the 
final technical report of Time- Zero, Inc. , on this nroject. 
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2. COSMIC- RAY NEUTRON MONITOR FOR OGO-VI 

The objective of this experiment is to monitor the 
cosmic-ray neutron flux over a large region of space near the 
earth for an extended time with a simple neutron monitor. The 
basic detecting system to measure the integrated neutron flux 
from 10“ 8 to 15 rieV consists of a moderated He 3 proportional 
counter surrounded by a charged-particle rejection system. 

A plastic scintillator viewed by a photomultiplier is operated 
in coincidence with the neutron counter to determine the 
energy spectrum and flux in the 1-10 MeV range. In order to 
minimize contributions from neutrons produced within the 
OGO spacecraft, the sensor is located about 4.5 meters out 
on the EP-5 boom. The experiment includes an in-flight 
calibrator to check pulse discriminator levels on neutron and 
charged-particle channels and a calibration- loop system to 
compensate for changes in the photomultiplier gain. The mean 
neutron detection efficiency for the exoected neutron energy 
spectrum over the energy range 10“ 8-15 MeV is 1.8 cm^. The 
weight of the neutron sensor is 2.58 kg, with a length of 
0.314 meter and a diameter of 0.111 meters. The total power 
requirements are 4.5 watts. 

(IEEE Transactions on Geoscience Electronics, 

Volume GE-7 , Number 2, April 1968, pp. 88-93.) 

The details of the design and testing of the detector are 
contained in Appendix 1, page A-l. 


3. LATITUDE AND ALTITUDE DEPENDENCE OF THE COSMIC 
RAY NEUTRON LEAKAGE FLUX 


Preliminary measurements of the cosmic ray albedo neutron 
flux above 400 km are reported for the period June 7 to 17, 
1969. The measurements were made with a detector on board the 
OGO-6 satellite that responds primarily to neutrons in the 
energy range 10 4 to 10 6 ev. The latitude variation of the 
counting rate was found to be 7.4/1 between 90° and 0°, 
similar to that predicted by Lingenfelter (1963). The 
neutron counting rate near the poles was found to decrease 
by 26+5% between 450 and 950 km altitude, corresponding to a 
r-4.2 ± 0.9 dependence, where R is the distance from the center 
of the earth. This large altitude dependence excludes any 
angular distributions of neutron flux at the top of the 
atmosphere that are more peaked toward the vertical than cos 
6 (where 9 is the angle from the vertical), but it is not in 



disagreement with an isotropic neutron flux. The total neutron 
leakage flux found by assuming isotropy was about 0.7 times 
that predicted by Lingenfelter (1963). The neutron leakage 
flux values obtained from the QGO 6 experiment agree with 
those estimated from other experiments with a similar neutron 
energy response, if the same angular distribution is used. 
Estimates of the total neutron leakage flux, deduced from 
experiments responsive only to neutrons above 1 MeV by 
using the Lingenfelter (1963) and Newkirk (1963) energy 
spectra of leakage neutrons , show agreement with the lower 
energy results when the Newkirk spectrum is used but are up 
to twice as great as the lower energy results when the 
Lingenfelter spectrum is applied. 

(Journal of Geophysical Research, Volume 75, No. 22, 
page 4197, 1970.) 

The details of the results on the latitude and altitude 
dependence of the cosmic ray albedo neutron flux are 
contained in Appendix II, page A2. 


4. ENERGY DEPENDENCE OF THE COSMIC RAY NEUTRON 
LEAKAGE FLUX IN THE RANGE .01 TO 10 MeV 

The cosmic-ray neutron leakage flux and energy spectrum 

in the range 1-10 iieV were measured by a neutron detector 

on the OGO 6 satellite from June 7 to September 30, 1969. 

The same detector simultaneously measured the total leakage 

flux, having 75% of its response to the leakage flux in the 

interval 1 keV to 1 MeV. For a neutron energy spectrum of 

the form AE Y i n the range 1-10 MeV, the upper limit to y for 

polar regions (P c <0.3 G v ) was found to be 1.0 and for the 

equatorial regions (P > 12 Gv) was 1.2. For the polar regions, 

the lower limit to y v)as found to be 0.8. This energy spectrum 

at 1-10 MeV is slightly flatter than L. L. Newkirk predicted. 

The neutron fluxes at 1-10 MeV were 0.28+0.03 and 0.035+0.003 

cm -2 sec -1 for P <0.3 and P > 12 Gv, respectively. The 

ratios of the nSutron flux c at 1-10 MeV to the total neutron 

flux were 0.54+0.07 for P > 0.3 Gv and 0.50+0.06 for P >12 Gv. 

c c 

(Journal of Geophysical Research, Vol. 76, Mo. 31., 

Page 7470, 1971.) 

The details of the analysis of the energy dependence of the 
neutron leakage flux are contained in Appendix III, page A3. 
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5. .UPPER LIMIT TO THE 1 TO 20 MEV SOLAR NEUTRON FLUX 

The upper limit on the quiet time solar neutron flux from 
1-20 MeV has been measured to be less than 2x10“ ^ n/cm^sec at the 
95% confidence level. This result is deduced from the OGO-VI 
neutron detector measurements of the "day-night" effect near the 
equator at low altitudes for the period from June 7, 1969 to 
December 23, 1969. The OGO-VI detector had very low (<4%) count- 
ing rate contributions from locally produced neutrons in the 
detecting system and the spacecraft and from charged-particle 
interactions in the neutron sensor. 

(Submitted to Solar Physics for publication) . 

The details of these searches for neutrons are contained in 

Appendix IV, page A4. 

6. A SEARCH FOR SOLAR NEUTRONS DURING SOLAR FLARE EVENTS 
During the operational period Q f the experiment from 
June 7 to December 23, 1969, several large solar flares 
occurred. No solar neutron flux was positively identified 
in the time association with large flare events in this 
time period. Hence, only upper limits to the solar neutron 
flux can be discussed. The results have been compared with 
the models for the solar neutron production during the solar 
flares proposed by Lingenfelter (1963) and the Lingenfelter 
and Ramaty (1967) . The details of the observations of the 
optical activity for the flares considered are in Table 5 and 
the solar proton emission in Table 6. Proton data were avail- 
able for all the solar flares. -For the large event on November 2, 
the data from Masley (1972) enabled us to determine the form 

of the energy spectrum using the exponential rigidity 

4 2 

function. J Q was 8.50xl0 4 protons/cm sec ster. , and 

2 

(dJ/dP) D = J Q /P o = 1850 protons per cm sec ster. M/v. 

P Q was approximately 50 MV. 



TABLE 5 


OPTICAL ACTIVITY 


Optical 
Flare (1969) 

Impor- 

tance 

Position 

Onset 
Time (ut) 

End 

(ut) 

Time of 
Max(ut) 

June 13 

3B 

S24E69 

15.49 

20.00 

16.33 

June 15 

2N-2B 

SI 717 7 7 

08.31 

08.50 

08.40 

Sept 25 

2N-3N 

N14W14 

07.00 

08.53 

07.24 

Nov 02 

3B 

N22W90 

10.28 

11.57 

11.39 

Dec 19 

in- 2 r 

N10W08 

12.01 

12.40 

12.20 


TABLE 6 


SOLAR PROTON EMISSION (EXPLORER 41) 



Solar 

Proton 

Flux 

Peak 

Date (1969) 

(ster 

cm 2 sec) 

-1 

time 


>60MeV 

>30MeV 

>10MeV 

Sept 25 

0.2 

1 

15 

-vll.OO 

Nov 02 

39 

257 

1437 

12.00-13.00 

Dec 20 

0.6 

1.1 

8.3 

'vOl.OO 
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The data used for this study were the gated neutron count 
rates for the OGO-vI detector on the sunward side of the earth. 

No background corrections were made because only the differences 
in counting rates were used. Only neutron rates corresponding 
to charged-particle rates less than 1.5 times the mean for that 
vertical cutoff rigidity were included in the analysis. To 
exclude solar-proton produced leakage neutrons we restricted 
the solar neutron analysis to geomagnetic latitudes less than 
30°. The neutron data were also selected from altitudes less 
than 900 kilometers to reduce neutron production effect by the 
energetic protons trapped in the inner radiation belt. 

Lingenfelter and Ramaty (1967) proposed that solar neutrons 
would arrive at the earth monoenergetic in time if the neutrons 
were emitted impulsively from the sun in a time interval which 
is comparable with the rise time of the optical flare. For these 
assumptions they calculated the time dependence neutron flux at 
the earth considering that neutrons will be produced both during 
the acceleration and the slowing down phases of the solar flare. 
The data which we selected were on the basis of the time depen- 
dency of the solar neutron flux as predicted by this model. 

For the large solar flare event on November 2, we found no 
statistically significant difference in the counting rates during 
the flare event and another pass through the same region when 
there would be no solar neutrons. For this event the upper limit 
to the solar neutron rate is 5.0x10 counts/sec at the 95% con- 
fidence level. The corresponding upper limit to the solar 
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2 

neutron flux is 0-135 neutron/cm per sec. If these results are 

corrected for the atmospheric solar neutron leakage flux arising 

from the interaction of solar neutrons with the earth's atmos- 

-2 2 

phere the upper limit is now 5x10 neutrons/cm per sec in the 
1-20 MeV energy range. 

Comparing this upper limit with the predicted solar neutron 

flux using the Lingenfelter model (1969), we find that the 

calculations predict a neutron flux at earth at approximately 

2 

0.21 neutrons/cm per sec. Noting that only about 27% of all 

the solar neutrons predicted by Lingenfelter are in the 1-20 MeV 

energy range, the maximum major upper limit to the total solar 

2 

neutron flux from our measurements is 0.19 neutrons/cm per sec. 
These observations, therefore, are consistent with this model. 

In other words , this upper limit is within agreement with 
Lingenfelter ' s prediction. This does not imply that there is 
a detectable solar neutron flux. For the other solar flare 
events on September 25, December 19, June 13, and June 15, 
upper limits to the solar neutron flux can also be obtained, 
not inconsistent with theoretical predictions. This we believe 
is the first time that an upper limit to a solar neutron flux 
has been measured for a flare of importance greater than 3B. 

We must emphasize that our measurements do not suggest that there 
are any solar neutrons from impulsive flare events. Perhaps 
only when larger detectors with better efficiency and direction- 
ality are flown will it be possible to determine if there are 
any solar neutrons. 



- 8 - 


7. SOLAR PROTON NEUTRON LEAKAGE MEASUREMENTS 

Energetic solar protons incident On the earth's atmosphere 
at high latitudes will produce neutrons, some of which will 
produce neutrons, which will then leak out of the atmosphere. 

A fraction of these leakage neutrons will travel in directions 
such that the radioactive decay of the neutrons into protons 
and electrons occurs within the magnetosphere. Consequently, 
the protons and electrons can become trapped. Hence, this may 
be one source of the energetic protons in the radiation belts. 
The neutron data from the OGO-VI satellite enabled us to 
evaluate the solar proton production of neutrons during several 
flare events . 

The neutron counting rates were derived from those events 
in the 3 He proportional counter not associated with the 
charged-particle guard counter events. These gated neutron 
counting rates were selected for geomagnetic latitudes X 
(north and south) leakage (greater than 70°) because the solar 
protons giving rise to the neutron are restricted to the polar 
regions of the earth. In most flare events the protons have 
kinetic energies less than about 300 MeV. Consequently, the 
effects of the solar protons will not be seen at latitudes less 
than 70°. The neutron counter rates were then corrected for 
dead time and reduced to a common altitude of 450 km using the 
observed altitude dependence discussed in Section 3. The 
details of the optical observations, riometer absorptions, and 
the solar proton emissions are summarized in Tables 1 and 2 
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for the solar flares of September 25, November 24, December 18, 
and December 19. Coincident with the proton flux enhancements 
recorded on Explorer 41 are the transient increases in the 
gated neutron rates above the cosmic ray albedo neutron rate 
background (Figures 1, 2, 3, 4) . 

To compare our results with the calculations of Lingenfelter 
and Flamm (1964) on the neutron leakage flux produced by solar 
protons we must have data on the energy spectra and charge 
composition of the energetic solar particles. Furthermore, we 
must make background corrections to the observed neutron 
counting rates. This also necessitates data on the solar par- 
ticle composition and energy spectra. The method to make 
these background corrections is straight-forward and is dis- 
cussed by Jenkins et.al. (1971). As an example of the background 
corrections we tabulate the contributions to the neutron rates 
of the 3 He counter at the poles for the September 25, 1969, 
solar flare event (Table 3) . 

After making the necessary background corrections for 
local production and non-neutron events, we obtain 0.1810.05 
counts/sec as the actual increase of the neutron counting 
rate at 450 km for the solar event on September 25. This 
corresponds to a neutron counting rate of 0.2310.06 counts/sec 
at the top of the atmosphere (taken to 50 km altitude) . There- 
fore, the neutron leakage flux from solar protons is 0.1110.03 
2 

neutrons/cm sec. From the measurements of the solar proton 
energy spectrum (Mas ley, 1972) we obtained a relation between 
the integral proton flux J and the rigidity of the protons P 



TABLE 3 

CORRECTIONS TO THE He 3 -COUNTER GATED 
NEUTRON COUNTING RATE (CTS/SEC) 

AT THE POLES FOR THE SEPTEMBER 25, 1969 
SOLAR FLARE EVENT 


Locally produced neutrons 

In the spacecraft mainbody 0.0052 

In the EP5 set up 0.0016 

In the neighboring electronics and 

detector walls 0.0576 

In the moderator 0.0009 

Highly ionizing events in the He 3 counter 

Low-energy solar protons 0.0000 

Low-energy solar a particles 0.0000 

Z >2 solar cosmic rays (including showers) 0.0000 

Evaporation protons from solar cosmic ray 

interactions 0.0004 

Proton and a-initiated shower particles 0.0000 


Total Correction 


0.0657 
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which agreed with the Explorer 41 proton data using the relation- 
ship J = J Q exp (-P/P Q ) and (dJ/dP) Q . We found that the 
characteristic rigidity of the solar protons, P , was 57 Mv 

and the total directional intensity of solar protons near the 

2 -1 

earth J was 60.4 (cm sec ster) . For this solar proton 
o 

rigidity spectrum the calculations of Lingenfelter and Flamm 

2 

(1964) predict a neutron leakage flux of = 0.085 neutrons/cm 
sec. Our result, therefore, is in reasonable agreement with 
the theoretical predictions for this flare event. 

We also studied the November 24, December 18, and December 19, 
solar proton events. The results are listed on Table 4. For 
these flares, however, the energy spectrum and charge composition 
of the solar particles were not available. We, therefore, 
assumed a power-law energy dependence and used the proton data 
of Explorer 41 on the integral flux greater than 60 MeV and 
greater than 30 MeV to estimate the integral and differential 
energy spectra. From Table 4, we see that the measured neutron 
leakage fluxes for the November 24, December 18, and 
December 19 solar-proton events are in reasonable agreement with 
the theoretical predictions. 

The measurement of the solar- proton neutron leakage at 
low latitudes is important because only the neutron decay 
products injected at low latitudes can contribute to the inner 
zone protons. The data used for the low latitude study were 
the gated neutron rates at geomagnetic latitudes both north 
and south, less than 45°. This selection minimizes the con- 
tributions from local production of neutrons by charged particle 
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interacting with the satellite assembly and unshielded portions 
of the detector. The neutron counting rates corresponding to 
charged particle rates greater than 1.5 times the minimum value 
for the particular cutoff rigidity were excluded from the 
analysis. The rates were selected for small intervals of 
latitude and altitude to avoid any latitude or altitude bias 
for the time intervals involved. There were no measureable 
Fbrbush decreases or other time-variations in the cosmic-ray 
source which would produce uncertainties in the interpretation 
of the results. As an example we show in Figure 5 the neutron 
rates for the September 25, 1969 event. There was no apparent 
enhancement in the neutron counting rates at low latitudes during 
the event. On the other hand for November 2, 1969 (Figure 6) 
there appears to be an enhancement in the neutron rates which 
is time associated with the solar-proton flux at high latitudes. 
Following the analysis of Lingenfelter and Flamm (1964) we can 
deduce the predicted neutron flux at low latitudes. The 
predicted neutron leakage rate at 42° geomagnetic latitude and 
750 km altitude is 0.21 counts /sec. This result is within the 
two standard deviations of the observed rate. E*r all the other 
events the predicted solar-proton neutron leakage fluxes at low 
latitudes were within the statistical uncertainties of the data, 
which for the known solar proton flux and energy spectra, agree 
with the theoretical predictions. 

Now let us investigate the possibility that the brief 
neutron increase on November 2 is from sources other than the 
solar proton albedo neutrons. First, the results on the 
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solar neutron flux (Section 5) indicate that it cannot be 
due to this source. No detectable solar neutron flux has 
been observed. Second, local production could not have 
produced the neutron enhancement because the charged particle 
rates were essentially the same as had been observed during 
previous passes through the same region. Third, any cutoff 
rigidity change from magnetic field perturbations can not be 
a cause for the enhanced neutron rates. Finally, it is un- 
likely that alpha particles could have made any significant 
neutron contributions . This is based on known proton-to-alpha 
ratios for the November 2 solar flare event. Therefore, for 
the November 2 event we must conclude that the brief neutron 
enhancement was presumably due to polar cap protons producing 
neutrons which then travelled to lower latitudes. 
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8. THE ANGULAR DISTRIBUTION OF THE NEUTRON LEAKAGE FLUX 
AS DEDUCED FROM THE ALTITUDE DEPENDENCE 


The altitude variation of the neutron counting rate above 
the atmosphere is related to the angular distribution of the 
leakage neutrons as follows: The counting rate of the neutron 

detector above the earth's atmosphere neglecting neutron decay 
is given by: 

f 

(1) N = ld(S) e (5) A (ft) 

s. ° 

4 it 2 

where e (ft) is the efficiency of the detector (counts/neutron/cm ) 

for neutrons arriving in the direction r at the detector. The 

neutron flux 4> (ft) in the direction ft at the spacecraft is the 

same as the neutron flux <}> (ft) at the top of the atmosphere if 

the direction (-ft measured from the detector) intersects the top 

of the atmosphere. Otherwise the flux is zero. We can write 

the neutron counter rate as follows 


(2) 


= / e 


N =/ dA e (ft) <t> (9) k (v) 


cosv 


where the integration is over the top of the atmosphere visible 
to the detector and <J> has been written as a separable function 
in 0 and v. These relationships are illustrated in Figure 7 
where we have indicated the spacecraft location and a sample 
point at the top of the atmosphere. In evaluating the above 
equation the variables were written in terms of height h above 
the top of the atmosphere, the geomagnetic colatitude and geomag- 
netic longitude. (The detailed calculations are presented in 
the Ph.D. dissertation by S. Ifedili) . In this analysis we must 
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keep in mind that the latitude dependence of the counting rate 
is much larger than the altitude dependence of the counting rate 
at any given latitude. Therefore the effects of latitude mixing 
upon the altitude dependence are very important. 

From the mean latitude and altitude dependence of the 
counting rates we can deduce the form of the angular distribu- 
tion for the neutron leakage in the energy range from .01 - 10 Mev. 
The neutron counting rates were again derived from the 3 He 
proportional counter events not associated with counts in the 
charged-particle guard counters. Contributions from local pro- 
duction were minimized by taking neutron counting rates corres- 
ponding to charged-particle rates less than twice the minimum 
values for the particular vertical cutoff rigidity. Data 
acquired in the regions near the Cape Town and Brazilian 
anomalies were excluded as well as data from excursions into 
the horns of the outer radiations belt. 

We must recognize that any error in the rigidity location 
or in the vertical cutoff rigidity might contribute an error to 
the altitude variation of the neutron intensity. Hence, we 
evaluated such contributions carefully. However, if each location 
in a rigidity bin contributed equally at all altitudes the errors 
would be a minimum. This is generally the situation when the period 
of analysis equals the time taken for the satellite's orbit to 
progress through 360°. This was approximately 4 months for the 
OGO-VI satellite. So the time period selected for the study 
was June 7 to September 30, 1969. The results of the altitude 
dependence are summarized in Table 7. The changes in the altitude 
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dependence of the neution leakage flux with latitude are attributed 
to the shape of the latitudb variation df the neutron flux and 
the effects of latitude mixing, recognizing that the detector 
sees neutrons from a wider range of latitudes as the altitude 
increases. We must remember that the neutron leakage flux 
starts to increase rapidly with latitude (Appendix II) at about 
20 °, increases steadily until about 55° when it flattens and is 
latitude independent above 55°. Hence, we would expect that at 
about 20° latitude, the higher altitudes, which include a larger 
range of geomagnetic latitudes, would be affected by the latitude 
mixing. At about 20° we would expect a minimum altitude change 
in the neutron counting rate. Similarly we expect a maximum 
altitude change for counting rates at locations greater than 55° 
when the counting rates are independent latitudes. This is 
exactly as observed. 

The percentage change in the counting rate as a function 
of geomagnetic latitude is shown in Fig. 8. The calculated 
altitude changes are indicated for different assumed angular 
distributions of the neutron leakage. The results shown suggest 
that the angular distribution of the neutron leakage at the top 
of the atmosphere is essentially isotropic. The best fit of the 
observed to the predicted altitude dependence at A <60° is for 
an isotropic angular distribution k(v)=l. Angular distributions 
which are peaked more to the vertical or to the horizontal do not 
agree as well with the observed altitude dependence. The better 
statistics in the polar regions enable us to make a better com- 
parison of the calculated and measured altitude dependence of the 
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neutron rates. For A >60° the angular distribution is given 
by k(v)= 1 -0.5 cos v, although an isotropic angular distri- 
bution. is not in disagreement with the observations. 

The somewhat flatter (or pancake-shaped) angular distribution 
at the poles might be due to the enrichment of the neutron 
energy spectrum preferentially in the high energy region. These 
high-energy neutrons would be very strongly peaked in the forward 
direction, being largely the result of knock-on collisions in 
the thicker atmosphere for which the neutrons travel near the 
horizontal direction. However, they have to be degraded in 
energy to less than 10 Mev in order to be detected by the 
OGO-VI sensor. From a survey of elastic high-energy neutron 
interactions, considerable peaking of the scattered neutrons 
in the forward direction is observed. Therefore, mixture of 
an isotropic high-energy neutrons degraded to less than 10 Mev 
(about 4%) with a preponderant fraction of evaporation neutrons 
(about 96%) might explain the flatter angular distribution in 
the polar region. 

9. SOLAR MODULATION OF THE COSMIC RAY NEUTRON LEAKAGE FLUX. 

The solar modulation of the cosmic rays should be reflected 
in the production of neutrons in the atmosphere of the earth. 

The effects of such a modulation would be observed in the rapid 
transient decreases (Forbush decreases) and in the 11-year 
variation which is anti correlated with the 11-year sunspot cycle. 
For a typical cosmic ray detector such as a neutron monitor the 
amplitude of the 11-year variation is about 20-30% and is energy 
dependent. The low energy particles are more strongly modulated. 
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Lingenfelter (1963) calculated the energy, latitude and 
altitude distributions for the neutron leakage at the minimum 
and maximum of the solar activity cycle. 

To check these predictions for the 11-year variation we 
must make certain that the data being used are free from all 
other variations. There should be no measureable lorbush 
decreases, no polar-cap neutron events, no observable solar 
neutrons, and no neutron production affects from radiation-belt 
protons. For the study we selected the months of July and 
October, 1969. During this period the Mt. Washington neutron 
monitor increased about 4%, which is one- fifth of the complete 
solar cycle variations. The OGO-VI data used were the 3 He gated 
neutron counting rates at 400-500 km altitude, excluding data 
acquired in the regions of the Capetown and Brazilian anomolies. 
The neutron counting rates were used only when the associated 
total charged particle rates were less than twice the minimum 
at that location. The neutron rates were also corrected for 
dead time, for locally produced neutrons, and highly ionized 
charged-particles using the techniques described in Appendix II. 

The results are shown in Figure 9. Assuming that any 
latitude the change in the neutron leakage flux as calculated 
by Lingenfelter was a linear function of the cosmic-ray intensity 
as measured by a neutron monitor at low cutoff rigidities, we 
could predict the neutron leakage rates for the months of July 
and October, 1969. We normalized the predicted neutron leakage 
flux to the observed 3 He count rates for July. From Figure 9 
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it is evident that the measurements are in reasonable agreement 
with Lingenfelter ' s calculations for the 11-year solar modulation, 
although the change for this time period is small. It appears 
that Lingenfelter ' s calculations can be used with much greater 
confidence to correct the observed neutron fluxes for the 11-year 
modulation which is necessary for comparison of the neutron 
fluxes measured at different epochs in the solar-activity cycle. 

We attempted to determine the effects of Forbush decreases on 
the cosmic-ray neutron leakage. Again the same criteria were 
applied in selecting the 3 He events. But, as Forbush decreases 
generally last only a few days, it is very difficult to make 
accurate comparisons with the large statistical errors present 
in the neutron data. A series of two successive Forbush 
decreases which occurred in November, 1969, provided an oppor- 
tunity to make limited comparisons. From this comparison we 
concluded that the cosmic-ray neutron flux (E<10Mev) at X>60° 
was depressed slightly during the Forbush decreases. However, 
the short duration of the Forbush decrease indicated that inter- 
grated effects on the cosmic-ray neutron leakage are insignifi- 


cant. 
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11. FIGURE CAPTIONS 


o 

- Neutron rates at geomagnetic latitude X>70 
during the September 25 solar proton event. 

o 

- Neutron rates at X>70 during the November 24 
solar nroton event. 

- Neutron rates at X>70° during the December 18 
solar nroton event. 

- Neutron rates at X>70° during the December 19 
solar nroton event. 

- Neutron rates at low geomagnetic latitudes 
during the September 25 solar nroton event. 

- Neutron rates at lovz geomagnetic latitudes 
during the November 2 solar nroton event. 

- The geometry for calculating the angular deoen- 
dence of the neutron leakage flux from the 
altitude variation of the observed neutron 
counting rate. 

- The angular dependence of the cosmic-ray 
neutron leakage flux. 

- The 11-year modulation of the cosmic-ray 
neutron leakage flux. 
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tion of phoswiching and pulse-shape discrimination 
(PSD) to separate neutrons from charged particles and 
gamma rays. The neutron energy spectrum itself was 
not determined and the results must be extrapolated to 
the top of the atmosphere to evaluate the neutron leak- 
age flux. The high-energy spectrum in the range 12 <£„ 
<100 MeV at atmospheric depths of 9 g-cm _! and 
geomagnetic latitude 42°N has been measured by Zych 
and Frye [4]. They observe that the proton recoil spec- 
trum was found to depend upon E ~ 3 - 3 . This proton 
spectrum is interpreted as due to a highly anisotropic, 
high-energy neutron flux near the top of the atmosphere. 
The extrapolated flux value at 0 g-cm 2 is 0.11 »/cm ! 
seconds, an order of magnitude greater than Lingen- 
felter [5] calculates. The most precise fast-neutron 
measurements (3<£»<20 MeV) have been made re- 
cently by St. Onge [6] with a proton recoil scintillator 
utilizing a separate charged-particle rejection scheme 
with a multiparameter display system for the PSD 
which separates gamma rays from neutrons. Unfor- 
tunately, data, so far, have been obtained only at the 
Pfotzer maximum. These results indicate that the neu- 
tron spectrum depends upon E~ a , where a is ~4 at E 
= 3.0anda~1.5 at 15 MeV. Also, a much larger ratio of 
gamma-rays to neutrons was observed, which suggests 
that the experimental difficulties are greater than antic- 
ipated. The neutron leakage flux in the energy interval 
0-15 MeV has been measured recently by rocket flights 
of moderated He 3 neutron detectors at several latitudes 
[7]. For latitudes of 0° and 70° the integral fluxes were 
0.08+0.02 and 0.80 + 0.10 «/cm 2 -s, respectively. These 
results are in better agreement with the neutron leakage 
fluxes calculated by Lingenfelter [5] and Newkirk [8] 
than earlier experimental values, although the latitude 
dependence of the integral flux is not so large as pre- 
dicted by Lingenfelter [5 ] . 

From this brief survey of the current knowledge about 
the neutron flux and energy spectrum both within and 
outside the atmosphere, we can see that many unsolved 
problems exist. First, within the atmosphere (at ~5 
g-cm -2 ) the neutron energy spectrum in the range 1 
<£»<50 MeV has not been precisely determined. Very 
few attempts have been made to determine the direc- 
tionality of the fast-neutron flux. Second, above the 
atmosphere there is disagreement between theory and 
experiment about the magnitude of the total neutron 
intensity. The energy spectrum in the 1-50 MeV range 
has not been precisely determined. The neutron in- 
tensity dependence upon distance from the earth is 
known only approximately [9] and definitive measure- 
ments have not been made of the angular distribution of 
the neutrons diffusing from the atmosphere. Both the 
theoretical and experimental values for the neutron 
leakage flux have been used to evaluate the contribu- 
tions which are made to the trapped radiation from the 
decay of these neutrons [10 ] , [l 1 ]. The results of such 
calculations indicate: 1) solar neutrons and solar pro- 
tons make an insignificant contribution; 2) the flux of 


trapped protons with £<20 MeV cannot be accounted 
for by the neutron decay mechanism; and 3) the magni- 
tude of the contribution from this source to the trapped 
protons with £>20 MeV is unresolved, the results from 
the theoretical calculations disagree with those expected 
from the measured leakage flux by more than an order of 
magnitude. Measurements of the fast neutron flux and 
energy spectrum should help resolve these differences. 

The experiment which we have designed can deter- 
mine the latitude dependence of the integral neutron 
flux in the energy range 10 -8 -15 MeV, as well as the 
longitude and altitude dependence. However, since the 
neutron sensor has an omni-directional response, no in- 
formation on the anisotropy of the neutron flux can be 
acquired. The neutron energy spectrum in the interval 
1-10 MeV can be deduced from the plastic scintillator 
used as a proton recoil detector in coincidence with the 
He 3 counter. 

Neutron Sensor 

A. Description 

Since the neutron intensity is measured in a back- 
ground of energetic charged particles, which can pro- 
duce neutrons in the sensor, it must include a scheme to 
reject neutrons associated with charged-particle events. 
Also, the sensor must clearly separate neutrons from 
gamma rays. The best estimates indicate that the inte- 
grated neutron leakage fluxes at latitudes of 90° and 0° 
are <1.0 and <0.08 »/ cm 2 -s, respectively. Conse- 
quently, the sensor must have a reasonable efficiency 
for neutron detection and a low background counting 
rate. 

The neutron sensor chosen for this mission is similar 
to the types flown on three USAF' Discoverer satellites 
[12] and, recently, on a series of rocket flights [7 ]. A He 3 
proportional counter, which detects neutrons indirectly 
through the nuclear reaction He 3 (n, p)H 3 , is surrounded 
by a plastic scintillator as shown in Fig. 1. The plastic 
scintillator acts both as a moderator for fast neutrons 
and a detector of proton and carbon recoils produced in 
the slowing down of the neutrons. This represents a 
unique application of a scintillator-type neutron de- 
tector. Hence, this moderated He 3 counter can be used 
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tern using laboratory beams of energy 2.4 and 14 MeV. 
In the actual flight unit the same definition of the pulse- 
height distribution cannot be obtained because only 
four pulse-height channels are used. These four channels 
are determined by the integral discrimination levels A- 
E inclusive. Crude measurements indicate that the 
efficiencies for the four channels (A-B, B-C, C-D, D- 
E) with the integral discriminators set as in Fig. 2 are 
0.03, 0.08, 0.14, and 0.05 cm 2 , respectively, at 2.4 MeV 
and are 0.013, 0.011, 0.017 and 0.030 cm 2 , respectively, 
at 14 MeV. 

In order to use the data in these four channels to de- 
termine the unknown energy spectrum the efficiencies of 
these channels for neutrons of many different energies 
from 1 to 20 MeV must be determined. The counting 
rate for each channel may be written 

Cj = £ Pj(Ei)N(Ei)AiE, for j = 1, 2, 3, 4 (1) 

i=i 

where j designates the channel number, Pj(Ei) = the 
efficiency of channel j for neutrons of energy E { , and 
A,E=energy interval, and N(E t ) = the unknown neu- 
tron energy spectrum. This assumes that the efficiencies 
Pj(Ei) have been determined for a large number of 
different neutron energies so that over A,E the values of 
Pj(Ei) and N(E) are essentially constant. If we assume 
that the energy spectrum is of the form N(E)=AE~r 
then (1) becomes 

C, = A £ Pj(Ei ) (£,"’) A,2Z for; = 1, 2, 3, 4. (2) 

i—i 

We know Cj from the experimental data and Pj(Ej) and 
A ( E from the calibrations. In order to find A and y we 
can assume values for them and calculate the expected 
values of Cj. A comparison between the calculated and 
measured Cj will determine what values, if any, of the 
assumed neutron spectrum parameters, A and y, are 
reasonable fits to the data. Other functional forms for 
the neutron spectrum may also be attempted in order to 
determine which form is best capable of fitting the data. 

C. Summary of Sensor Properties 

The specifications for the neutron sensor are as fol- 
lows: weight = 5.7 lb (2.58 kg); length = 12f inches. 
(0.314 meter); diameter = 4f inches (0.111 meter); 
power requirements 4.5 watts. The estimated efficiency 
for the neutron leakage spectra is about 1.8 cm 2 . In the 
fast mode the efficiency at any energy value in the in- 
terval 1-10 MeV is about 10 percent of the efficiency in 
the total mode. In addition to the leakage flux measure- 
ments, the energy spectrum 1-10 MeV can be evaluated 
from the other simultaneous mode of operation. The 
sensor has excellent discrimination against gamma rays 
-and a low background counting rate. 


D. Electronics System 

In discussing the electronics system for the experi- 
ment, we will consider first that section pertaining to the 
charged-particle monitoring, then the section of the sys- 
tem related to the total neutron flux, and finally those 
parts dealing with the fast-neutron measurements. 
Since the experiment required more data output than 
could be handled with only the three main frame words 
provided us on the spacecraft, an internal subcommuta- 
tor sequence was incorporated within the experiment. 
In Fig. 3, only the experiment words are indicated 
with no reference to the subcommutator sequence. 

Referring to Fig. 3, we see that both coincidence 
rates A-C and B-D and the total counting rate (A+B 
+ C + D) of the charged-particle guard counters is read 
out on separate experiment words designated words 2, 3, 
and 5, respectively. Since extremely high counting rates 
(>10 4 per second) are expected when the spacecraft 
passes near and through the radiation belts, a logarith- 
mic counting rate meter is provided (spacecraft sub- 
commutator analog word 86). Since it is expected that 
the charged -particle rates, both total and coincidence (as 
well as the total neutron rate) will vary from 1 to 10 4 
per second, a logarithmic prescaling system was incor- 
porated into these channels (experiment words 2, 3, and 
5). If the counting rate is low so that less than 64 events 
are stored between readouts, then the number of events 
accumulated equals the word reading. But if, for ex- 
ample, the word is greater than 63 but less than 128, the 
number of events equals 65 +2 X (word reading 64). 
Thus, the maximum information which can be stored is 
then 2756 for experiment word 5 (9-bit word) and 708 
for experiment words 1, 2, and 3 (8-bit words). 

The neutron counter output is directly fed into experi- 
ment word 1, designated the ungated or total counting 
rate. This rate is due to neutrons incident from outside 



EXP'T WORD 2 
EXP'T WORD 3 
ANALOGUE WORD (SC 86) 
j-EXP'T WORD 5 

EXP'T WORD 4 
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h-EXP T WORD 6 


ANALOGUE WORD 
(SC 54) 


Fig. 3. Schematic diagram of electronics system. 
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Preliminary measurements of the cosmic ray albedo neutron flux above 400 km are reported 
for the period June 7 to 17, 1969. The measurements were made with a detector on board the 
Ogo 6 satellite that responds primarily to neutrons in the energy range 10* to 10’ ev. The 
latitude variation of the counting rate was found to be 7.4/1 between 90° and 0°, similar to 
that predicted by Lingenfelter (1983). The neutron counting rate near the poles was found 
to decrease by 26 :£ 8% between 450 and 950 km altitude, corresponding to a /£-*•» de- 
pendence, where R is the distance from the center of the earth. This large altitude dependence 
excludes any angular distributions of neutron flux at the top of the atmosphere that are more 
peaked toward the vertical than cos $ (where 6 is the angle from the vertical), but it is not in 
disagreement with an isotropic neutron flux. The total neutron leakage flux found by assum- 
ing isotropy was about 0.7 times that predicted by Lingenfelter (1963). The neutron leakage 
flux values obtained from the Ogo 6 experiment agree with those estimated from other 
experiments with a similar neutron energy response, if the same angular distribution is 
used. Estimates of the total neutron leakage flux, deduced from experiments responsive only 
to neutrons above 1 Mev by using the Lingenfelter (1963) and Newkirk (1963) energy spectra 
of leakage neutrons, show agreement with the lower energy results when the Newkirk spec- 
trum is used but are up to twice as great as the lower energy results when the Lingenfelter 
spectrum is applied. 


In this paper we present preliminary meas- 
urements of the cosmic ray albedo neutron flux 
in the energy range lO’-lO' ev from a neutron 
detector on board the orbiting geophysical ob- 
servatory Ogo 6, a polar-orbiting satellite with 
altitudes between 400 and 1100 km. From these 
measurements, the latitude and altitude depend- 
ence of the neutron flax for the initial turn-on 
period, June 7-17, 1969, have been determined. 
The June 7-17 period is typical in terms of the 
solar and geomagnetic activity present during 
this phase of the current solar cycle, with no 
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exceptional solar flare events reported [ Solar 
Geophysical Data, 1969]. 

Detector Description 

The components of the neutron sensor and its 
orientation with respect to the spacecraft and 
the earth are depicted in Figure 1. The sensor 
was located approximately 5 meters out on the 
EP-5 boom to minimise the background of 
neutrons that are produced in the main body 
of the spacecraft by the cosmic radiation. 

J The detector consists of a proportional counter 
filled with He”, sensitive primarily to thermal 
neutrons through the reaction He*(n, p)H”, 


4197 



4198 


JENKINS, LOCKWOOD, IFEDILI. AND CHUPP 




a-c a s-o 



Fig. 1. Neutron sensor arrangement and orientation on Ogo 6 spacecraft. 


surrounded by a plastic scintillator and a 
ring of charged-particle guard counters. The 
plastic scintillator serves two functions. First, 
it acts as a moderator for fast neutrons, by 
slowing them down to energies where they can 
be detected more efficiently by the He" counter. 
Second, it is used to measure the neutron energy 
in the range 1-10 Mev. The neutron energy 
distribution is deduced from the pulse height 
distribution of scintillation events followed 
within 25 microseconds by a He* counter event. 
The ring of charged-particle guard counters acts 
as an anti-coincidence shield against cosmic 
ray nuclei that can produce neutrons in the 
material adjacent to the He* counter. As shown 
in Figure 1, the charged-particle counters are 
arranged in four banks, A, B, C, and D. 
The signals from the four banks are summed 
and are also put into the coincidence pairs A-C 
and B-D. A more complete description of the 
sensors with associated electronics has been 
published [ Lockwood et al., 1969]. 

The following data are acquired by the neu- 
tron sensor: (a) the total number of events in 
tiie He* proportional counter (N T ) ; (6) the 
number of events in the He* counter not asso- 
ciated with events in the charged-particle guard 
counters (N a ) ; (c) the number of scintillator 
ervents associated with events in category (b), 
stored in a 4-channel pulse-height analyzer; 


(d) the number of events in the charged-par- 
ticle guard counter (N AB CD ) ; (e) the numbers 
of coincidences between the two pairs of 
charged-particle counter banks (N M , N„ B ). In 
this paper we are reporting total neutron meas- 
urements and concern ourselves primarily with 
the quantities Nr, N a , N ASC d, N iC , and N tB . 

Detector Efficiency 

The neutron detector was calibrated as a 
function of neutron energy by exposing it to 
known neutron fluxes from both Am“‘-Li, 
Am*°-Be, and Am^-moek fission radioactive 
sources and monoenergetic neutrons ranging 
from 5 kev to 20 Mev produced in the reactions 
Sc"(p, n) Ti w , H*(p, n)He*, H*(d, n)He“, and 
H*(d, n)He\ by using the 5.5-Mev Van der 
Graaf accelerator at the. Oak Ridge National 
Laboratories. The neutron detector efficiency 
for events in the He* counter, as a function of 
energy for an isotropic neutron flux, is shown 
in Figure 2. 

The measurements of detector efficiency for 
neutrons from the radioactive sources provided 
a check on the absolute efficiencies deduced 
from the monoenergetic measurements. From 
the observed variation of the monoenergetic 
neutron efficiency with energy and the neutron 
energy distributions for the Am“-Be, Am“-Li, 
and Am*“-mock fission neutron sources [Han- 
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sen, 1960; Geiger and Hargrove, 1964] , the 
efficiencies expected for these radioactive 
sources were calculated. The measured efficien- 
cies for the sources were then plotted in Fig- 
ure 2, at neutron energies for which the mono- 
energetic curve gives the calculated efficiencies. 
It can be seen from Figure 2 that the efficiencies 
measured for the radioactive sources agree with 
those obtained from the monoenergetic neutron 
measurements. 

To extrapolate the neutron efficiency for 
the detector to lower energies, we have used 
the measurements of Tatruta et d. [1965] with 
cylindrical paraffin-moderated BF» counters. 
The solid curve in Figure 2 is the expected 
efficiency as a function of energy based on 
their measurements. To adapt the results of 
Tatsuta’s measurements to the Ogo 6 detector, 
we interpolated his results to derive a neutron 
efficiency versus energy curve for an amount of 
hydrogen per unit length of moderator equal 
to that of the Ogo 6 detector. The neutron 
efficiency curve thus derived was normalized to 
the measured values obtained in the energy 
range 5 kev to 2.5 Mev. It is evident from 
Figure 2 that the deduced curve agrees with 
the energy dependence of efficiency found from 
0.1 to 2,5 Mev. At energies above 2.5 Mev, the 
efficiencies measured here are greater than 
those deduced from Tatsuta’s measurements, 
probably due to the moderating effects of in- 
elastic neutron collisions with the relatively 
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greater number of carbon nuclei surrounding 
the Ogo 6 detector. 

To derive a mean efficiency for the albedo 
neutron flux, the efficiency curve (Figure 2) 
was folded into the albedo neutron energy 
spectrum predicted for the present level of 
solar activity by Lingenjelter [1963]. The 
mean efficiency was also derived by using the 
calculated neutron energy spectrum of New- 
kirk [1963], The mean efficiencies obtained for 
an isotropic neutron flux and the anticipated 
contributions from the various neutron energy 
ranges to the counting rate of our detector are 
listed in Table 1. From this table we see that 
the largest fraction of the counting rate i8 from 
neutrons in the energy range 10* to 10“ ev for 
both energy spectra. The mean efficiencies for 
the two assumed spectra of leakage neutrons 
differ by 22%. 

The variation of the detector efficiency with 
the incident angle of the neutron flux was de- 
termined by using neutrons from the Am*”-Li 
source {(En) = 390 kev) and the reaction 
H“(d, n)He* (En = 14 Mev). The angular de- 
pendence of relative response shown in Figure 
3 is energy independent. The measured angular 
response curve (Figure 3) was used to convert 
the measurements taken with monoenergetic 
neutrons incident perpendicular to the axis of 
symmetry of the detector into the absolute 
efficiencies for isotropically incident neutron 
fluxes plotted in Figure 2. Because the re- 
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TABLE 1. Predicted Contributions to the Count- 
ing Rate of the Ogo 6 Neutron Detector* 


Neutron Percentage of Neutron Counting Rate 
Energy, Lingenfeller Newkirk 

ev [1963] Spectrum [1963] Spectrum 


<10 l 

2.3 

4.7 

10*-10> 

4.2 

6.0 

10M0* 

6.7 

8.0 

10M0* 

10.8 

13.4 

lOMO* 

16.3 

15.2 

10M0» 

47.4 

33.8 

10M0 1 

11.5 

16.3 

>W 

0.8 

2.6 

Mean detector 



efficiency for 



isotropic flux, 



cm* 

1.59 ± 0.16 

1.30 ± 0.13 


* From various energy ranges of leakage neutrons 
for the Lingenfeller [1963] and the Newkirk [1963] 
energy spectra for neutron leakage flux. 


sponse of the detector depends upon the direc- 
tion of the incident neutrons, the neutron flux 
in space derived from the observed counting 
rate will depend on the assumed angular dis- 
tribution of neutrons in the vicinity of the 
detector. The effect is small for the orientation 
of the neutron detector on the Ogo 6 space- 
craft, however, because all physically accept- 
able angular distributions are cylindrically sym- 
metrical about the vertical. The calculated 
response for different angular distributions 
varies by less than 15% from that expected 
for an isotropic neutron flux. 

Correction for Locally Produced Neutrons 

The total neutron counting rate N T of the 
Ogo 6 detector arises from three sources: (1) 
neutrons external to the spacecraft, i.e., not 
produced in the spacecraft or detector; (2) neu- 
trons produced in the detector by charged 
particles seen by the guard counters; and (3) 
neutrons produced in the spacecraft and de- 
tector by charged particles not seen in the 
guard counter. The gated neutron counting rate, 
N 0 , arises from neutrons in categories (1) and 
(3). The counting rate due to the neutrons ex- 
ternal to the experiment ( N a ') was estimated 
by u&ng the increase in neutron counting rate 
when the Ogo 6 satellite passed through the 
South Atlantic magnetic anomaly. Protons 
with energies of the order of 100 Mev that were 


present in the anomaly gave rise to enhanced 
local neutron production. From these observa- 
tions the ratio of counts from locally produced 
neutrons in category (3) to counts from cate- 
gory (2) neutrons was obtained. During other 
times, the gated neutron rate N 0 could be cor- 
rected for neutrons in category (3) by noting 
the count rate due to neutrons in category (2), 
(N r — No), and by using the ratio found in 
the anomaly. We found that the counting rates 
due to locally produced neutrons in categories 
(2) and (3) were approximately 35% and 12% 
of the total neutron counting rate N r , respec- 
tively. 

Results 

The gated neutron data corrected for locally 
produced neutrons (N 0 ') were sorted according 
to geomagnetic latitude and altitude into 5® by 
100-km bins. Neutron counting rates corre- 
sponding to charged particle rates in excess of 
2.5 times the minimum values for that latitude 
were excluded from this analysis. 

Latitude and altitude dependence. In Fig- 
ure 4 the neutron rates N 0 ' for the 400- to 
500-km altitude region are plotted as a function 
of geomagnetic latitude. The latitude depen- 
dence of counting rate from the poles to the 
equator is found to be 7.4/1, in good agreement 
with that predicted by Lingenfelter [1963] for 
this period in the solar activity cycle. The 
mean efficiency of the Ogo 6 neutron detector 
was assumed to be independent of latitude, al- 
though there was a 2% variation introduced 
from changes in the energy spectrum of leakage 
neutrons with latitude. 

The change in neutron rate N 0 ' with the al- 
titude can be expected to depend upon the 
latitude at which the measurements were taken. 
This is a result of latitude mixing, i.e., neutrons 
from a wider range of latitudes contribute to 
the counting rate as the altitude is increased. 
Since the altitude variation as seen by the Ogo 
6 detector is much smaller than the latitude 
variation, the latitude variations must be cor- 
rectly evaluated before determining the altitude 
variation. The data presented here used the 
predicted latitudes for the Ogo 6 spacecraft, 
which may be in error by several degrees. Such 
systematic errors can greatly affect the altitude 
variation of neutron rate N„ at latitudes where 
the neutron flux changes rapidly with latitude. 
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Fig. 3. Angular dependence of efficiency of detector. 


For these reasons we limit our present examina- 
tion of the altitude dependence of neutrons to 
latitudes greater than 70“ , where the neutron 
flux can be seen to be constant over latitude 
(Figure 4). 

The altitude dependence of Ay for these lati- 
tudes is shown in Figure 5. The mean decrease 
in count rate Ay between 450 and 950 km is 
26.0 ±: 4.6% (rms error). Written in terms of 
a radial dependence of neutron count rate, this 
becomes N a '{R ) = AR~ l -‘ * where R is the 
distance from the center of the earth. For com- 
parison, an altitude variation in Ay expected 
for a completely vertical neutron flux (R'“ ° al- 
titude dependence) is shown in Figure 5. 
A vertical neutron flux disagrees with our 
observations (probability of chi-squared greater 
than observed is 0.005). Our observations are 



Fig. 4. Latitude dependence of neutron count- 
ing rate Ns' for altitudes between 400 and 500 
km. 


not in disagreement with the altitude variation 
R' 3 ° between 450 and 950 km altitude expected 
for an isotropic neutron flux at the top of the 
atmosphere (probability of chi-equared greater 
than observed is 0.16). Taking the significance 
level in chi-squared to be 0.05, our observations 
agree with angular distributions of neutron flux 
that are flatter than cos 6, where 6 is the angle 
from the vertical. For our purposes here, we 
have assumed an isotropic neutron flux at the 
top of the atmosphere. 

Neutron leakage flux ai the top of the at- 
mosphere. To determine the leakage flux at 
the top of the atmosphere, the observed neu- 
tron rates between 400 and 500 km, denoted 
AV ( 450 ), were reduced to the top of the at- 
mosphere (taken to be 50 km) by using the 
calculated altitude dependence expected for an 
isotropic neutron flux distribution at the top 
of the atmosphere. This amounted to a mean 
correction of 43%. Then the counting rate at 
50 km, Ay (50), was converted to the leakage 
flux fa, by using the relationship fa = 0.5 Ay 
(50) /(e) for a 2ir-isotropic neutron flux distri- 
bution, where («> is the mean efficiency. This 
same relationship was used by Haymes [1964] 
and Holt,et al., [1966], who had detectors re- 
sponding primarily to 1- to 10,-Mev neutrons. 
Other experiments using detectors sensitive pri- 
marily to neutrons of lower energy [ Bame 
et al., 1963; Boella et al., 1965; Lockwood and 
Friling, 1968] have determined the neutron 
flux <j> — N/{(). This is equal to the leakage 
flux only for the case of an angular distribution 
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TABLE 2. Total Leakage Neutron Fluxes Found from the Ogo 6 Experiment and Several 1-to 10-Mev 

Experiments. 


Latitude, 


Flux Found with 

Flux Found with 

degrees 

Experiment 

Lingenfelter Spectrum* 

Newkirk Spectrum* 

7.4 

Ogo 6 

0.06 ± 0.01 

0.07 ± 0.01 

8 

Holt et at. [1966] 

0.09 

0.06 

42.3 

Ogo 6 

0.25 ± 0.03 

0.30 ± 0.04 

41.5 

Haymes [1964] 

0.36 ± 0.03 

0.24 ± 0.03 

42 

Holt et at. [1966] 

0.36 

0.24 

68.1 

Ogo 6 

0.42 ± 0.05 

0.52 ± 0.06 

68 

Holt el at. [1966] 

0.84 

0.64 


* The Lingenfelter [1963] and Newkirk [1963] energy spectra for leakage neutrons. 


in which all the neutrons are moving vertically 
upward at the top of the atmosphere. As seen 
in the preceding section, such an angular dis- 
tribution is not compatible with the observed 
altitude dependence in neutron counting rate. 

In Table 2 are given total neutron fluxes for 
the Ogo 6 detector, computed from the rates 
No' by using the mean efficiencies for the Ling- 
enfetter [1963] and Newkirk [1963] energy 
spectra of leakage neutrons. Also shown are 
total neutron fluxes obtained from measure- 
ments in the 1- to 10-Mev region by Haymes 
[1964] and Holt ct at., 1968]. The 1- to 10-Mev 
measurements have had several corrections ap- 
plied to them, to compare them with the Ogo 6 
results. First, the values were corrected to the 
same period in the solar-activity cycle by using 
the Mt. Washington neutron monitor rate and 
the calculated relative leakage fluxes between 
solar ma ximum and minimum [Ling enfetter, 
1963]. Second, the 1- to 10-Mev neutron fluxes 
were multiplied by the ratio of total to 1- to 
10-Mev neutrons given by the appropriate neu- 
tron energy spectrum (Lingenfelter or New- 
kirk) to determine a total neutron leakage flux. 



Fig. 6. Altitude dependence of neutron counting 
rate No for latitudes > 70°. 


It can be seen from Table 2 that the total neu- 
tron fluxes found from the 1- to 10-Mev meas- 
urements and the lower energy Ogo 6 measure- 
ments are in closer agreement if the Newkirk 
spectrum of leakage neutrons in the vicinity of 
spectrum is used. If the Lingenfelter spectrum 
is used, the 1-to 10-Mev measurements yield 
fluxes that are as much as twice the correspond- 
ing Ogo 6 flux, thus indicating that the energy 
1 Mev is harder than that calculated by Lingen- 
felter. 

In Figure 6 the values of total neutron leak- 
age flux <j> L , found for the Ogo 6 detector by 
using the Newkirk energy spectrum of leak- 
age neutrons, are plotted as a function of 
geomagnetic latitude. 

The other experimental values (Figure 6) 
were corrected for solar modulation to com- 
pare them with the present results. They were 
also corrected for altitude by using the al- 
titude dependence expected for an isotropic 
neutron flux distribution at the top of the 
atmosphere. Where necessary, they were ad- 
justed to the neutron leakage flux dependence 
on counting rate given by 4>l — 0.5 N/{t). Fin- 
ally, the Newkirk energy spectrum has been 
used with the published efficiency versus energy 
curves for these experiments to derive the 
neutron flux from the counting rates. 

The total neutron leakage flux measurements 
reported for the Ogo 6 detector are in agree- 
ment with the other experimental observations 
(Figure 6). Also, the measurements reported for 
the 1- to 10-Mev experiments agree with the 
lower -energy measurements when the Newkirk 
energy, spectrum is used. This is not true for 
the Lingenfelter energy spectrum (Table 2). 
Spectral measurements in the 1-to 10-Mev 
region [Haymes, 1964; Holt et at., 1966] have 




COSMIC RAY ALBEDO NEUTRON FLU'X 4203 



Fig. 6. Neutron albedo flux as a function of latitude. 


found an energy spectrum similar to or slightly 
harder than the Newkirk spectrum. 

The total neutron leakage flux as a function 
of latitude given by Lingenfelter [1963] plotted 
in Figure 6 gives the same calculated intensity 
as that of Newkirk at 57°, where the Newkirk 
calculation was made. The experimental obser- 
vations of total leakage flux in Figure 6 are 0.7 
times the predicted values. This result is 
consistent with the observations of Miles [1964] 
on the neutron densities in the atmosphere, 
which gave values 0.6 times that derived from 
the Lingenfelter calculations. 

Summary 

The present results regarding the neutron 
leakage flux indicate the following : a total 
neutron leakage flux of 0.7 times the intensity 
predicted by Lingenfelter [1963] and Newkirk 
[1963] ; a latitude variation similar to that 
given by Lingenfelter [1963] ; an energy spec- 
trum below 1(P ev similar to that given by 
Newkirk [1963]; an angular distribution of 
neutrons at the top of 'the atmosphere which 
may be isotropic over 2v but in any case is 
unlikely to be more peaked toward the vertical 
than a cos 6 distribution, where 6 is the angle 
from the vertical. 
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The cosmic-ray neutron leakage flux and energy spectrum in the range 1-10 Mev were measured 
by a neutron detector on the Ogo 6 satellite from June 7 to September 30, 1969. The same detector 
simultaneously measured the total leakage flux, having 75% of its response to the leakage flux 
in the interval 1 kev to 1 Mev. For a neutron energy spectrum of the form AE~ y in the range 
1-10 Mev, the upper limit to y for polar regions ( P c < 0.3 Gv) was found to be 1.0 and for the 
equatorial regions (P c > 12 Gv) was 1.2. For the polar regions, the lower limit to y was found 
to be 0.8. This energy spectrum at 1-10 Mev is slightly flatter than L. L. Newkirk predicted. 
The neutron fluxes at 1-10 Mev were 0.28 ± 0.03 and 0.035 ± 0.003 cm - * sec -1 for P c ? 0.3 
and P. > 12 Gv, respectively. The ratios of the neutron flux at 1-10 Mev to the total neutron 
flux were 0.54 ± 0.07 for P„ > 0.3 Gv and 0.50 ± 0.06 for P, > 12 Gv. 


The energy spectrum of the cosmic-ray neu- 
tron leakage flux in the range 1-10 Mev has 
,been measured by several workers [Haymes, 
1964; Baird and Wilson, 1966; Holt et al., 
1966; Albemhe and Talon, 1969; and Merker 
et al., 1970] who used detectors sensitive to 
neutrons through the production of knock-on 
protons. These measurements have been com- 
pared with lower-energy measurements made at 
other times [Jenkins et al., 1970] in order to 
extend the knowledge of the energy spectrum 
to the range 0.001-10 Mev. In the present 
paper we report the results of four months of 
simultaneous measurements of leakage neutrons 
in the ranges 1-10 Mev and 0.001-1 Mev made 
with the University of New Hampshire neutron 
detector on the polar-orbiting Ogo 6 satellite. 
From these data, information is deduced about 
the energy spectrum in the range 0.001-10 Mev 
and its variation with latitude. 

Detector 

S The neutron detector used has been described 
> in detail by Lockwood et al. [1969] and Jenkins 
et al. [1970] , It consists of a “He-filled propor- 
' tional counter, which 13 sensitive mainly -to slow 
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neutrons through the reaction ’He (n, p)*H 
and is surrounded by a plastic scintillator and 
a ring of charged-particle guard counters. The 
plastic scintillator serves two functions; first, 
to slow down faster neutrons so that they can 
be detected more efficiently by the *He counter; 
and, second, to measure the neutron energy 
spectrum in the range 1-10 Mev with a four- 
channel pulse-height analyzer. The neutron 
energy distribution is deduced from the pulse- 
height distribution of scintillation events fol- 
lowed within 25 jisec by a 'He-counter event. 
The charged-particle guard counters act as an 
anticoincidence shield against cosmic-ray nuclei 
that may produce detectable neutrons in the 
material adjacent to the ’He counter. Data 
obtained include the total counting rate of the 
“He-filled counter (TN), the rate of “He-counter 
events not gated off by the charged-particle 
guard counters (GN), and the rates of events' 
in the four scintillation channels that precede 
GN events by les3 than 25 /usee (PHA 1, 2, 3, 
and 4). 

The efficiencies of the “He counter (GN) and 
four scintillator channels (PHA 1, 2, 3, and 4 ) 
were measured at energies from 5 kev to 20 
Mev with a 5.5-Mev Van der Graaf accelerator. 
The efficiency is plotted in Figure I as a func- 
tion bf energy for the “He-counter events. (The 
details of the efficiency calibrations are given by 
Jenkins et al. [1970].) Energy spectra calcu- 
lated for the cosmic-ray leakage neutron flux 
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Fig. 1. “He neutron counter efficiency as a function of energy. The details on these 
measurements are presented by Jenkins et al. [1970]. 


1. Polar region, where the cosmic-ray ver- 
tical cutoff rigidity P. at the earth’s 
surface is <03 Gv and altitude is be- 
tween 400 and 500 km. 

2. Equatorial region, where P. is >12 Gv 
and altitude is between 400 and 600 km. 

The counting rates GN and PHA 1, 2, 3, and 
4 were corrected for a small background arising 
from neutrons produced locally in the space- 
craft and surrounding detector material by 
cosmic rays, and also arising from highly ioniz- 
ing cosmic rays whose energy loss in the “He 
counter is above threshold. The present back- 
ground correction represents a revision of the 
correction applied in a previous paper [ Jenkins 
et d., 1970], However, both the present cor- 
rection. (5% for GN) and previous correction 
(12% for GN) are small enough. that the con- 
clusions of that paper are not changed. The 
details of the background correction are given 
'in Appendix 1. 

- After corrections for background, the rates 
GN and PHA 1, 2, 3, and 4 were reduced to the 
top of the atmosphere (taken to be at 50-km 
altitude), using the altitude dependence calcu- 


[Newkirk, 1963; Lingenfelter, 1963] were used 
to deduce mean efficiencies for the “He counter 
[Jenkins et d., 1970]. From these efficiencies 
tiie observed counting rates (GN) can be con- 
verted to a total neutron leakage flux. For 
either the Newkirk or Lingenfelter spectra, the 
“He counter was found to respond mainly to 
neutrpns between 10 kev and 1 Mev. 

The efficiency as a function of energy for the 
scintillator channels PHA 1, 2, 3, and 4 is 
shown in Figure 2. Although the upper dis- 
crimination levels were set below 10-Mev energy 
loss, there is a finite efficiency for neutrons of 
energies above 10 Mev because of the contribu- 
tion of inelastic carbon collisions to the slowing 
down of higher-energy neutrons in the scintil- 
lator. Uncertainties in the measured efficiencies 
of PHA 1, 2, 3, and 4 are 18%, 16%, 8%, and 
Tl%, respectively. 

Data Reduction 

The data were obtained between June 7 and 
September 30, 1969, and have been sorted ac- 
cording to spacecraft position. In this paper we 
are concerned principally with two regions: 
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lated for our detector assuming a 2 ir isotropic 
angular distribution of leakage neutrons at tlie 
top of the atmosphere. From the counting rates 
at the top of the atmosphere and the measured 
efficiencies, we then calculated the fluxes of 
leakage neutrons, again assuming an isotropic 
angular distribution for the neutrons leaking 
out of the atmosphere. 

Spectral Shape in ’rap Range 1-10 Mev 

The scintillator rates PHA 1, 2, 3, and 4 were 
used to examine the shape of the neutron leak- 
age spectrum at 1-10 Mev. An energy spectrum 
of the form AE~ y was assumed and was folded 
together with the scintillator efficiency curves 
of Figure 2 to derive the expected counting 
rates of PHA 1, 2, 3, and 4 for various values 
of A and y. The expected and observed rates 
were then compared in both the polar and 
equatorial regions by a x' test to determine 
which values of A and y best fitted the obser- 
vations. The results of this test for the polar 
region data in June 1969 are shown in Figure 3. 
Contour lines representing boundaries to the 
regions within which 60% and 90% of the 



Fig. 2. Neutron efficiencies as a function of 
energy for the four scintillator channels (PHA 1, 
2, 3, and 4). 



Fig. 3. Fit of AE^r to the neutron PHA results 
for June 1969 at P, < 0.3 Gv and altitudes <500 
km. All data reduced to 50 km, using an altitude 
correction to the neutron counting rates. Contour 
line for boundaries to regions within whioh 60% 
and 90% of allowed values must fall are shown. 

allowed values of A and y must fall are shown. 
From this graph it is clear that any spectral 
shape flatter than E~ l fits our data. This is a 
consequence of the finite efficiencies of four 
scintillator channels at neutron energies above 
10 Mev and the observed counting rates. The 
neutron leakage flux in the energy region 1-10 
Mev, computed from the allowed values of A 
and y for the polar region in June 1969, is 
0.28 ± 0.03 neutrons cm'* sec* 1 (Table 1). An 
upper limit (at the 95% confidence level) to 
“the steepness of the energy spectrum y m « = 
1.07 is found. The rest of the results for the 
four months of polar and equatorial regions 
are presented in Table 1. The neutron fluxes 
at 1-10 Mev in both regions for the four months 
agree within the given rms errors. Generally, a 


TABLE 1. The Upper limit to 7 (touui) and 1-10 Mev Neutron 
Leakage Flux (ATi.it * /i u AE’~'[ dB) 


Month 
(1969) . 

Polar Region 

Equatorial Region 

• Tarn* 

2Vi_te, cm” 1 sec* 1 

Tnmz 

Ni-it, cm"* see"» 

June 

1.07 

0.28 ± 0.03 

1.28 

0.032 ±0.003 

July ■ 

0.95 

0.28 

1.12 

0.035 

Aug. *' 

0.88 

0.28 

1.07 

0.036 

Sept. 

1.02 

0.27 

0.85 

0.032 
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TABLE 2. Comparison of Normalized Predicted Counting Rates and Observed Rates (sec -1 ) for the Four 
Scintillator Channels PHA 1, 2, 3, and 4 in June 1969 


Predictions Observations 


Channel 

Newkirk, 

1963 

Lingenjelter, 

1963 

Wilson et al., 
1969 

Merker, 

1970 


PHA 1 

0.0020 A: 0.0004 

Equatorial Region ( P c > 12 Gv) 
0.0026±0.0005 0.0019 ±0.0003 

0.0021 ±0.0004 

0.0020±0.0002 

PHA 2 

0.0026±0.0004 

0.0030 ±0.0005 

0.0025 ±0.0004 

0.0027 ±0.0004 

0.0022 ±0.0003 

PHA 3 

0.0022±0.0002 

0.0021 ±0.0002 

0.0023 ±0.0002 

0.0022 ±0.0002 

0.0019±0.0002 

PHA 4 

0.0029 ±0.0003 

0.0020 ±0.0002 

0.0030±0.0003 

0.0027 ±0.0003 

0.0034±0.0003 

X* (3 deg 
freedom) 

3.2 

16.8 

2.5 

5.1 


Significance 
of fit, % 

40 

<0.1 

45 

15 


PHA 1 

0.017±0.003 

Polar Region (Pc < 0.3 Gv) 
0.024±0.004 0.016±0.003 

0.019 ±0.003 

0.015±0.001 

PHA 2 

0.023±0.004 

0.028±0.004 

0.022±0.003 

0.024±0.004 

0.019±0.001 

PHA 3 

0.020 ±0.002 

0.020 ±0.002 

0.020±0.002 

0.020±0.002 

0.018±0.001 

PHA 4 

0.026±0.003 

0.019±0.002 

0.027 ±0.003 

0.024±0.003 

0.030±0.001 

X 1 (3 deg 
freedom) 

4.8 

34.6 

3.2 

9.0 


Significance 
of fit, % 

20 

<0.1 

35 

3 



lower y m „ is found in the polar than in the 
equatorial region, but this is not significant in 
light of the larger statistical errors in the equa- 
torial region. 

A second technique was used to test the 
published calculations of leakage neutron energy 
spectra. The spectral distributions calculated by 
Newkirk [1963], Lingenjelter [1963], Wilson 
et al. [1969], and Merker [1970] were folded 
together with the scintillator efficiencies and 
normalized so that the best fit to our observa- 
tions would be obtained. The normalized' pre- 


dictions and the observed rates (PHA 1, 2, 3, 
and 4) were then compared. The results for 
June 1969 are shown in Table 2. The Lingen- 
felter spectral shape, given by E~ UI> in the energy 
range 1-10 Mev, does not fit our observations. 
The Newkirk spectrum, which can be approxi- 
mated by E~ 10 in this energy interval, is a fair 
fit to the June data. The more recent calcula- 
tions of Wilson et al. [1969], in which some 
structure in the spectrum is evident, fit the 
June data possibly better than the Newkirk 
calculations. The Wilson calculation, however, 


TABLE 3. Fit of Newkirk [1963] and Wilson et al. [1969] Spectra to the Observed PHA Counting Rates for 
’ June-September 1969 


Month 

(1969) 


Equatorial Region 



Polar Region 


Newkirk 

Significance 

level* 

Wilson et al. 

Newkirk 

Significance 

level* 

Wilson et al. 

X* 

% 

x* , 

■ % 

x’ 

% 

X* 

% 

June 

3.2 

40 

2.5 • 

45 

4.8 

20 

3.2 

35 

July 

3.6 

! 30 

2.2 

55 

7.9 

5 

5.7 

15 

Aug. 

6.0 

10 

5.0 

. 20 

10.1 

1 

7.6 

5 

Sept. 

10.1 

1 

7.6 - 

5 

6.4 

10 

4.5 

20 


* 3 degrees of freedom. 
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TABLE 4. Comparison of the Fast Neutron Leakage Flux(ATi-io) with the Total Leakage ( N T ) using Newkirk'a [1963J Energy 
. Spectrum 





Polar Region 



Equatorial Region 




0889) 


N r 

N,/N r 

A7 

N r 


(N r /N r )«> 


Jane 

0.80 ±0.02 

0.56 ±0.06 

0.64 ±0.07 

0.034 ±0.003 

0.076 ± 9.0008 

0.46 ±0.06 

+0.090 ±0.025 


JTaly 

0.29 

0.55 

0.54 

0.039 

0.077 

0.51 

+0.028 


Aug. 

0.82 

0.66 

0.67 

0.040 

0.077 

0.53 

+0.040 


Sept. 

0.80 

0.68 

0.61 

0.039 

0.079 

0.60 

+0.012 


, : Predicted ratio at X 

~67°N 

0.46 






v eoittods only to 10 Mev. Contributions to the ratio of fast neutron leakage flux to total neu- 

; counting rate from above this energy were tron leakage , flux near the poles is 052, in 

ignored for this spectrum. Both these spectra contrast to the measured ratio for the June 

are slightly steeper than the observations would data of 0.68 ± 0.09. Similarly, for the equa- 

■: indicate. The Merker spectrum is a poor fit, torial region the calculated ratio is 0.32 and the 

being too steep for the observations. measured ratio (June) is 0.54 ± 0.07. Thus 

The fit of the observed PHA counting rates the Lingenfelter spectrum does not predict a 

to the Newkirk and Wilson et al. spectra is sufficiently high fraction of leakage neutrons 

generally poorer for the later months (Table 3). at 1-10 Mev for the present measurements. This 

The observed counting rates in the PHA chan- result supports the preliminary results [Jenkins 

nds for all four months agree within the et al., 1970] in which the 'He-counter measure- 

statistical fluctuations. Therefore we must con- ments were compared with the fast neutron 

(flude that the Newkirk and Wilson et al. spectra measurements made by other workers [Haymes, 

are only poor fits to the observed data. From 1964; Holt et d., 1966]. 

the relative counting rates in the four channels, Better agreement between the predicted and 
a flatter energy spectrum than those calcu- measured ratio of fast neutron leakage flux to 
lated is indicated for the later months as well. 


Comparison op Neutron Flux at 1-10 Mev 
with Total Neutron Flux 

Information about the energy spectrum of the 
leakage neutrons can be extended to lower 
energies by comparing the scintillator measure- 
ments (PHA 1, 2, 3, and 4) with the 'He-counter 
measurements (GN). From the normalization 
constants used in fitting the calculated spectra 
Of Newkirk [1983] and Lingenfelter [1963] to 
. the observed scintillator counting rates, leakage 
fluxes in the region 1-10 Mev (fast neutrons) 
obtained’ for each of these spectra. The 
4 V total leakage fluxes were obtained by convert- 
ing the *He-counter rate to a flux using the 
, corresponding neutron efficiencies deduced for 
each of the two spectra. The ratio of fast neu- 
tron leakage flux to total leakage flux thus 
found from , the scintillator and Tie-counter 
. rates was ; then compared with the ratios pre- 
• dieted directly from the Lingenfelter [1963] 
and Newkirk [1963] spectra. 

For the Lingenfelter spectrum, the predicted 



Fig. 4. -Differential neutron energy spectra. The 
Lingenfelter and Newkirk spectra have been 
normalized to the Ogo 6 measured neutron flux 
for E < L Mev. Above E = 10 Mev the spectrum 
is assumed to be of the form E~°. 
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Fig. 6. Geometry of the neutron detector used 
to calculate the neutron production in the de- 
tector. 


total neutron leakage flux is obtained using the 
Newkirk [1963] spectrum (Table 4). For the 
equatorial regions the Newkirk spectral shape 
agrees with our measurements. For the polar 
region the result suggests that there are rela- 
tively more fast neutrons than the Newkirk 
spectral shape (originally calculated for a geo- 
magnetic latitude of 57°N) predicts. The in- 
crease in the observed ratio of fast-to-total 
leakage flux from the equator to the poles is 
also listed in Table 4. A relatively small error 
arises because the large errors in detector effi- 
ciency contribute in the same way to errors 
in the ratio at both positions; i.e., the efficiency, 
error expressed as a fraction is constant. This 
enrichment of the fast neutron leakage spec- 
trum is seen in the calculations of Lingenfelter 
[1963] ; however, the observed change is 
greater. BoeUa et al. [1965] interpreted his 
observed variation of latitude effect with at- 
mospheric depth as being due to more of the 
. neutron sources being close to the top of the 
atmosphere at high latitudes, and' anticipated a 
resultant enrichment of the 1- to 10-Mev leak- , 
age flux at high latitudes. 

Finally, the measurements of the total neu- 
tron leakage flux and the neutron leakage flux 
at 1-10 Mev can be used in conjunction with 
information on the spectral shape for E < 1 
Mev and E > 10 Mev to check the measured 


values of y in the interval 1-10 Mev. For the 
polar regions, the minimum value of x* yields 
a 1- to 10-Mev leakage flux given by N(E) = 
0.06£k°“ neutrons cm'* sec' 1 . Let us consider 
how this spectrum fits the calculated spectra 
at lower energies. For the Newkirk neutron 
energy spectrum, 87% of the counting rate of 
the Ogo 6 *He detector is contributed by neu- 
trons with E < 1 Mev [Jenkins et al., 1970]. 
Therefore let us normalize the Newkirk spec- 
trum of E < 1 Mev by the ratio of the neutron 
flux measured by the Ogo 6 to the total flux 
at 0 g/cm* calculated by Newkirk. The New- 
kirk spectrum for 2? < 1 Mev normalized by 
this ratio (0.560/0.899) is shown in Figure 4. 
For the Lingenfelter spectrum, 81% of the 
counting rate of the- Ogo 6 detector arises from 
neutrons E < 1 Mev. Thus we can normalize 
the Lingenfelter neutron leakage spectrum for 
sunspot minimum and A. > 80° by the ratio 
of the measured total leakage flux to the calcu- 
lated flux (0.460/1.367) as shown in Figure 4. 
It is known that for E < 10 kev the spectrum 
goes as E" 1 , and the only differences in the two 
spectra are for E > 10 kev. However, at E = 
1 Mev, N(E) for the normalized spectra are 
equal. Consequently, the differential spectrum 
for 1 < E < 10 Mev is most probably given 
by 0.10 A' 0 80 , or as steep as 0.13 E~ 1M . Any 
spectrum as flat as E'°' a does not merge with 
either the Lingenfelter or Newkirk spectra for 
E < 1 Mev. For a flat spectrum the leakage 
flux for E > 10 Mev is much larger than that 
calculated by Lingenfelter. For either the nor- 
malized Lingenfelter or Newkirk spectra in 
Figure 4, the ratio of the leakage flux at 1-10 
Mev to the total leakage flux is 0.46, in ap- 
proximate agreement with the observed value. 
Therefore, for the polar regions the lower limit 


TABLE 5. Corrections to the 'He-Counter Gated Counting 
Rate (counts per second) et 750-km Altitude ior-P. dt 03 Gv 


Locally produced neutrons 

In the spacecraft , 0.003 

In the neighboring electronics and detector walls 0 .014 
In the moderator 0.003 

JEghly ionising events in the 'He counter 

Low-energy cosmio-ray protons 0.005 

'• Low-energy cosmic-ray a particles 0 .010 

Z > 8 cosmic rays (including showers) 0.008 

Evaporation protons from cosmic-ray interactions 0.001 
Proton and a-initlated shower particles 0.000 

Total correction (estimated uncertainty. 4-30%) 0.043 
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to y is 0.8, with the most probable value 
0.8 .< y < 1.0 for the interval 1-10 Mev. 

Summary 

We are able to set the limits on y in the 
polar neutron energy spectrum (AR' 1 , 1 < 

E < 10 Mev) as y m «. = 1.0 and y mlo = 0.80. 

This agrees with the results of Holt et al. 

[1966] and Merker et al. [1970], who found 
y = 1.05 ± 0.16 and y = 1.08 ± 0.10, respec- 
tively, near the top of the atmosphere in the 
same energy range. Our measured upper limit 
for y indicates a flatter spectrum than the value 
y = 1.3 ± 0.10 at 41.5°N measured by Haymes 
[1964]. The measurements by Baird and Wilson 
[1966] of the 1- to 10-Mev spectrum on a 
rocket yield y = 0.8 ± 0.3, essentially in 
agreement with our results. A slightly steeper 
spectrum (y = 1.24) was deduced by Albemhe 
and Talon [1969] for the range 3-14 Mev. 

However, no errors on y were given. There is 
fair agreement with the leakage spectrum calcu- 
lated by Newkirk [1963] and Wilson et al. 

[1969] in the fast-neutron range (1-10 Mev). 

The measured ratio of fast neutron leakage 
to' the total leakage increases slightly between 
the equator and poles. The observed ratio near 



Fig. 6. Background corrections to the counting rates of the “He counter (GN) and 
v the scintillator channels (PHA 1, 2, 3, and 4) as a function of vertical cutoff rigidity (P«) . 


the equator agrees with calculations by New- 
kirk [1963] done for a latitude of 57°N, but 
near the poles the ratio indicates that there is 
a greater fraction of fast neutrons that the 
Newkirk calculation predicts. 

Appendix: Corrections 

Background correction of the ‘H e-counter 
rate. The corrections to the “He gated count- 
ing rate that are needed before it can be used 
to determine the leakage neutron fluxes arise 
from the following types of events. 

a. Neutrons produced locally in the space- 
craft, adjacent electronic circuitry, and moder- 
ator by cosmic rays escaping detection in the 
anticoincidence guard counts. 

b. Charged particles entering the unguarded 
ends of the “He counter, which, by their high 
ionization loss in the “He counter, produce volt- 
age pulses above discriminator threshold. The 
following particle events fall into this classifica- 
tion: 

1. Cosmic-ray protons with original energy be- 
tween 100 and 125 Mev. 

2. Cosmic-ray a particles with original energy 
between 100 and 375 Mev/nucleon. 
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TABLE 6. Background Contributions (per cent) 
to the Counting Rate of the Scintillator PHA 
Channels for the Polar Region 


Particles 

PHA 1 

PHA 2 

PHA 3 PIIA 4 

Locally produced 

3.6 

3.8 

3.7 3.0 

neutrons 

Low-energy 

0.2 

0.5 

0.6 3.5 

protons 

Total 

3.8 

4.3 

4.3 6.5 


3. Cosmic-ray L, M, and H nuclei of all ener- 
gies sufficient to reach the “He counter. 

4. Evaporation protons produced in the moder- 
ator and walls of the “He counter by cosmic 
rays. 

5. Shower particles produced in nuclear inter- 
actions of the cosmic rays with the modera- 
tor and walls of the “He counter. 

Corrections to the “He-counter gated counting 
rate were evaluated for these sources with the 
cosmic-ray fluxes given by Webber [1967] cor- 
rected to 1 the present time in the solar activity 
cycle by using the Mount Washington neutron 
monitor intensity. The cosmic-ray cutoff rigidi- 
ties (P„) were evaluated from the tables given 
by Shea et al. [1968], The neutron production 
cross sections and resulting energy distributions 
of the locally produced neutrons were taken 
from Chen et al. [1955], Dostrovski et al. 
[1958], Jain et al. [1959], Jain [1961], Berco- 
vitch et al. [1960], Bertini [1965], and AUmiller 
et al. [1967], The geometry used for these cal- 
culations is shown in Figure 5. 

Partial corrections calculated for the polar 
region at 750 km altitude are listed in Table 5. 
The total correction as a function of geomag- 
netic cutoff is plotted in Figure 6. 

Background correction of the scintillator rates. 
A scintillator output pulse will occur for several 
of the types of “He-counter background events 
listed above. However, since the uppermost dis- 
crimination level on the scintillator corresponds 
to a 10-Mev proton energy loss, the geometry 
involved excludes a large fraction of the back- 
ground because the events correspond to more 
than a 10-Mev proton energy loss in the scin- 
tillator. Only a small fraction of the cosmic-ray 
protons counted by the “He detector as GN 
events have a sufficiently small scintillator light 
output to be counted as PHA events. The con- 


tribution from particles with Z > 2 is negligible. 
On the other hand, the locally produced neu- 
trons having energies between 1 and 15 Mev 
can contribute to the PHA events. 

The contribution to the counting rates PIIA 
1, 2, 3, and 4 of the scintillator from locally 
produced neutrons was determined by assuming 
an evaporation energy spectrum. The spectrum of 
Alsmiller et al. [1967] for evaporation neutrons 
produced by 400-Mev protons incident on 
aluminum was assumed. Then the relative contri- 
butions (JoN,fi , i = 1, 2, 3, 4) of these neutrons 
to the counting rates GN and PHA 1, 2, 3, and 
4 were evaluated from the measured efficiency 
curves. The absolute contribution of locally 
produced neutrons to the PHA 1, 2, 3, and 4 
rates was then determined by putting = 
(f </ S on) Ray, where Ray is the absolute contri- 
bution of background to the gated neutron rate. 

The low-energy proton contributions to the 
PIIA counting rates were deduced by noting the 
relative increases in GN and PHA 1, 2, 3, and 
4 rates that occurred in the South Atlantic 
anomaly, where largo fluxes of protons are 
known to be present. From these relative 
increases and the normal contribution of low- 
energy protons to GN calculated in the pre- 
ceding section, the normal contribution of low- 
energy protons to PHA 1, 2, 3, and 4 rates were 
evaluated. Table 6 lists the background correc- 
tions to the scintillator as percentages of the 
normal counting rate for the polar region. The 
corrections are less than 7% at all latitudes. 
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ABSTRACT 

The upper limit on the quiet time solar neutron flux 

- 3 

from 1-20 Mev has been measured to be less than 2 x 10 
n/cm sec at the 95% confidence level. This result is 
deduced from the OGO-6 neutron detector measurements of 
the "day-night" effect near the equator at low altitudes 
for the period from June 7, 1969, to December 23, 1969. 

The OGO«6 detector had very low(<4%) counting rate con- 
tributions from locally produced neutrons in the detecting 
system and the spacecraft and from ohar^ed-p^tAcLe 1q.1~.p-x*- 
actions in the oeutron sensor* 


1 


Now at Communications Research Center, 


Ottawa, 


Ontario, Canada 




1 . 


INTRODUCTION 


Biermann, Haxel, and Schluter (1951) first suggested 
that solar flare protons with energies of about 100 MeV 
could interact with solar atmosphere to produce a secondary 
flux of neutrons. Later theoretical calculations by Hess 
(1962), Lingenfelter et al . (1965a, b), and Lingenf elter 
and Ram at y (1967) gave further support to the possible 
existence of a solar neutron flux and showed that the only 
significant source of solar neutrons is from protons acc- 
elerated and slowed down during flares. For almost two 
decades efforts have been made to detect such a neutron flux. 
The measurement of a solar neutron flux is important because 
this is a means to evaluate some parameters of the accelerat- 
ing mechanism of solar cosmic— rays su-ch as the time scale 
for acceleration, the nuclear reactions which may occur, and 
the region in the solar envelope in which acceleration takes 
place. Attempts to measure this flux have been made both 
during small solar flares and during relatively quiet periods. 
The existence of a detectable solar neutron flux has not 
been confirmed in the searches by Hess and Kaifer <1967) , 
game and Ashbridge (1966), Webber and Ormes (1967), Forrest 
and Chyipp (1969), Zych and Frye (1969), Kim (1970), Heidbreder 
et al . (1970), Cortellessa et al. (1970, 1971), Daniel 


0.9 71), and Eyles et al . (1971). No attempt will be made 
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here to discuss all these results. For the details, 
reference should be made to the review of solar flare 
neutron production by chupp (1971) and the survey of the 
experimental efforts to detect solar neutrons compiled 
by Lockwood (1972). 

We report here the results of the 0G0-6 neutron 
experiment which set a new lower limit to the 1-20 MeV 
neutrons produced at the sun during a relatively quiet 
period. These new measurements will be compared with 
others in the energy region 1-100 MeV. 

2 . INSTRUMENTATION 

The neutron sensor used for these measurements is a 

3 

moderated He proportional counter surrounded by an 
anti-ooincideaco charged particle shield composed of 22 
proportional counters and has a solar neutron energy 
response from about 10 ^ MeV to 20 MeV. The moderator 
is a plastic scintillator which acts both as a moderator 
for fast neutrons and a detector of proton and carbon 
recoils produced in the slowing down of the neutrons. 

The neutron detecting system and its location on the 
OGO-6 spacecraft are shown in Figure 1. The details of 
the operation of the detector have been described else- 
where ( Lockwood et al ., 1969). 

The neutron detector was calibrated as a function of 
neutron energy by exposing it to known neutron fluxes 


Fig . 
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from both Am 2lfl -Li, Am 2l+1 -Be, and Am 2l+1 -mock fission radio 
active sources and monoenergetic neutrons ranging from 
5 keV to 20 MeV produced in reactions using the Van der 
Graaf accelerator at Oak Ridge National Laboratories. 

3 

The neutron detector efficiency for events in the He 
counter as a function of energy for an isotropic neutron 
flux is shown in Figure 2. The details of the detector 
calibration have been presented ( Jenkins et al . , 1970) . 

The variation of the detector efficiency with the inci- 
dent angle of the neutron flux was determined at neutron 
energies of 0.3 MeV and 14 MeV. This measured angular 
response curve was used to convert the measurements 
taken with monoenergetic neutrons incident perpendicular 
to the axis of symmetry of the detector into absolute 
efficiencies for isotropically incident neutron fluxes 
plotted in Figure 2. Since the response of the detector 
depends upon the direction of the incident neutrons, 
the neutron flux in space derived from ^he observed count- 
ing rate will depend upon the -assumed angular distribution 
of neutrons in the vicinity of the detector. 

3 . RESULTS 

The neutron counting rates were derived from those 
3 

events in the He proportional counter not associated with 
events in the charged-particle guard counters. This is 
referred to as the "gated" neutron counting rate. No 


FIG. 2 
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3 

correction was made for background events in the He 
counter because differences in counting rates rather than 

absolute counting rates were used. In any case this back- 

( 

ground rate was only about 3% at the equator ( Jenkins 
et al . , 19 ' 0 , 1971). Any contributions from locally pro- 
duced neutrjns in the 0G0-6 spacecraft were calculated to 
be less than. 0.3%. 

The data for June 7, 1969, to December 23, 1969, 

inclusive were first sorted according to spacecraft 

location. The minimum values of the charged particle 

counting rates for all vertical cutoff rigidities (P ) 

c 

and altitudes showed only a cutoff rigidity dependence. 

The ratio of the polar (P <0.5 Gv) to equatorial (P 

c c 

>14 Gv) minimum charged-particle rates was 4.3 to 1. "-> 

minimize the contributions from local charged-particle 
production the neutron rates corresponding to charged 
particle rates in excess of 1.5 times the minimum values 
for that rigidity were excluded in the analyses. The 
gated neutron counting rates were then sorted according 
to vertical cutoff rigidity and altitude into 18 vertical 
cutoff rigidity intervals for seven differ^:.'; altitude 
ranges. The calculated cutoff rigidity values 15 Shea et al . 
(1968) were used in this sorting. Data from the regions 
near the Capetown and Brazilian anomalies in t .e geomag- 
netic field were excluded from the analysis. we also 
eliminated data for those days on which olar particle 
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events were reported by detectors on Explorer 41, 

Pioneer 8 and Pioneer 9. The analysis for solar neu- 
trons was restricted to only the equatorial regions 
for which the vertical cutoff rigidity at the earth's 
surface is >14 Gv. T.'his restriction minimized the 
problems due to intensity - time variations and neutrons 
produced in the atmosphere at high latitudes by solar 
protons during any solar flare events which might not 
have been excluded. Finally, the neutron data were 
selected from altitudes less than 900 Km to remove neu- 
tron production effects by the energetic protons trapped 
in the inner radiation belt. 

The differences between the neutron counting rates 
for the detector on the sun-ward side ("day") and for 
the detector completely eclipsed by the earth ("night") 
are listed in Table 1 for low altitudes and high cut- 
off rigidities. The errors listed are statistical and 
are equal to the square root of the counting rate 
divided by the time of measurements. There ire no 
statistically significant differences between the two 
counting rates. From these data the weighted average 
solar neutron rate n g in the 1-20 MeV energy range is 


n = 


N 

2 

i=i 


n , 
l 


£ 


! N 

I v 1 _ 2 
a . 2 j 0 i 
-i + >/izi 


a . ‘ 

i 


TABLE I 


( 1 ) 



TABLE I, DIFFERENCES IN "DAY" AND "NIGHT" NEUTRON COUNTING RATES 


RIGIDITY 

AVERAGE GEO- 

ALTITUDE 




COUNTING 

RATES (/SEC) 


INTERVAL 
... (gv) 

HmG |~AT a 


(km) 


DAY 

NIGHT 

DIFFERENCE 

14 - is 

7.5 

400 

- 500 

0.14344 


0.00094 

0.14670 

± 0.00220 

- 0.0033 ± 0.0024 

> 16 

0 

400 

- 500 

0.1200 

4. 

0.0013 

0.1178 

± 0.0021 

+ 0.0022 ± 0.0024 

14 - is 

7 . 5 

500 

- 600 

0.1372 

+ 

0.0013 

0 . 1417 

± 0 . 0029 

- 0.0045 ± 0.0032 

> 16 

: • 0 

500 

- 600 

0 . 1129 

± 

0.0016 

0.1177 

± 0.0031 

- 0.0047 ± 0.0035 

14 - 16 

7.5 . 

600 

- 700 

• 0.1338 

i 

0 . 0014 

0.1340 

± 0.0025 

- 0,0002 ± 0.0029 

’> 16 

0 

600 

- 700 

0.1102 

+ 

0.0017 

0 . 1104 

± 0.0030 

- 0 . 0001 5 ± 0.0035 

14 - 15 

7.5 

700 

- 800 

0.1247 


0.0016 

0.1219 

± 0.0021 

+ 0.0028 ± 0.0026 

' > 16 

0 

700 

- 800 

0.1059 

+ 

0.0019 

0.1054 

± 0.0024 

+ 0.0005 ± 0,0031 

14 - 16 

7.5 

800 

- 900 

0.1197 


0.0018 

0.1176 

± 0.0019 

+ 0.0021 + 0,0025 

> 16 

0 

800 

- 900 

0.1012 

zh 

0.0022 

0.1003 

± 0.0022 

+ 0 . 0009 ± . 0 . 0030 
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= (-1.54 ± 9.30) x 10~ 4 /sec. 

Therefore, the upper limit to the continuous solar 

- 3 

neutron rate is l.£6 x 10 /sec at the 95% confidence 
level . 

The resulting upper limit continuous solar neutron 
flux F can be calculated from 



C is a factor which depends upon the assumed angular 

distribution for the neutrons incident on the detector. 

2 

The mean efficiency e in cm was obtained by folding 

the efficiency efE^) for neutrons incident perpendicular 

to the axis of symmetry of the detector into a solar 

neutron energy spectrum of the shape calculated by 

Lingenfelter and Ramaty (1967). The shape of the solar 

neutron energy spectrum at earth depends upon the solar 

charged-particle spectrum which is given by dJ/dP = 

(dJ/dP) Q exp(-P/P Q ) where P is the particle rigidity and 

P q is the characteristic rigidity which varies from event 

to event ( Freier and Webber , 1963). Differential solar 

neutron energy spectra for P Q ranging from 60 to 400 MV 

have been presented by Lingenfelter and Ramaty (1967). If 

this differential neutron spectrum is taken to be of the 

2 

form n(E) neutrons/cm sec MeV , then the mean efficiency 



-7 


is given by 

e = / n (E j e (E) dE//n (E) dE. (3) 

The mean efficiencies for spectra with P = 60 MV and 

o 

2 

P = 125 MV are calculated to be respectively 0.38 cm 
o 

2 

and 0.37 cm , about 20% larger than for an isotorpic flux. 

If neutrons are incident perpendicular to the axis of 

symmetry, then C = 1. Therefore, at the 95% confidence 

level, the upper limit to the integral solar neutron 

-3 2 

flux in the energy range 1-20 MeV is 4.9 x 10 /cm sec 

for P q = 60 MV. The neutron energy range is defined 

by the low-energy shape of the solar neutron energy 

spectrum which decreases rapidly at lower neutron 

energies because the survival probability of the low 

energy neutrons is small. 

If there are any solar neutrons, then this upper 

limit to the measured direct solar neutron flux must be 

corrected for the atmospheric neutron leakage flux 

arising from the interaction of solar neutrons with the 

earth’s atmosphere. This correction is carried out as 

follows. Let the total solar neutron count rate be 

N = N, + IJ , where the subscript A refers to the contri- 
A D 

bution from atmosphereic leakage and D to the contribution 
from solar neutrons directly incident on the detector. 
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The calculations of Alsmilief and Boughner (1968) for 
the differential neutron fixix (neUtrons/cm sec MeV) pro- 
duce at vafioUS depths in the atmosphere by solar flare 
neutrons can be used to estimate the resulting neutron 
leakage flux. Then, the corresponding solar neutron 
fluxes (see equation 2) in the energy interval covered 
by the OGO-6 detector are: F = N / Z and F, = N / z e, , 
since C = 2 for a flux isotropic over a hemisphere. 
Putting F = a F ' 


F 


D 




(4) 


Now, since the OGO-6 detector responds principally to 

direct solar neutrons with energies 1-20 MeV, F^ = F T 

(/ 20 J(E) dE//” J(E)dE) = 0.267F t . F t is the total 

integral incident solar neutron flux with a differential 

energy spectrum J(E) given by Lingenfelter and Ramaty 

(1968) for P q = 60 MV. The incident solar neutron 

flux F gives rise to an atmospheric neutron source F , 

T N 

at small atmospheric depths, of which a fraction 

T 

r = (F ) /F leaks out of the atmosphere. Of the total 

A JM 

T 

solar neutron leakage F ft , the OGO-6 detector responds 

to a flux F . from 10 keV-10 MeV. Putting b = F_/F_^, 

A A A 

then F (io keV-10 MeV) = b F /’ = b r F„ = brn F m = 

A ANT 

^ ^ ^ Fp = a F^. Assuming the neutron spectral shape 
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for the leakage flux to be the same as given by Hess 

e t al . (1961) , b = 0.94 at the geomagnetic equator for 

the period in the solar cycle corresponding to that 

for these measurements. For this spectrum r is estimated 

to be 5 < r < 10% ( Lingenf elter , 1963) . Taking the 

smallest value of n = 1 to yield the highest upper 

limit to F , a = 0.25. For this neutron leakage 

2 

spectrum = 1.35 cm . Therefore, the corrected upper 
limit to the quiet time solar neutron flux from equation 
(4) is 

F d = 0.363 = 1. 8xl0 _3 n/cm 2 sec. 

e D 

In evaluating the solar neutron leakage flux from 

the atmosphere the contributions from neutrons with 

<0.5 MeV were assumed to be given by the spectral 

shape of Hess e t al . (1961) . Actually, the spectrum of 

Alsmiller and Boughner (1968) is still riched in neutrons 

below 0.5 MeV. For the latter spectrum the upper limit 

is reached by about 10%. Such an extrapolation is not 

included in the upper limit placed on the solar neutron 

flux by the 0G0-6 measurements. 

The upper limit to the "quiet-time" integral solar 

neutron flux deduced here is comparable to the value of 
-3 2 

< 2 x 10 cm/sec in the 1-10 MeV range obtained by 
Hess and Kaifer (1967). Upper limits to the solar neutron 
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fitix have also been set by £ych arid Ftye (1969) , Kim 
( 1967 ), Forrest and Chupp (1968), and Cortellessa et al. 
(1971), but at much higher energies from 10 to 200 MeV . 

For comparison with other measurements the CG0-6 result 
was converted to a differential flux by folding the 
energy dependent response function into the shape of the 
solar neutron energy spectrum calculated by Lingenf el ter 
and Ramaty (1967) for a characteristic rigidity P q = 

60 MV. For characteristic rigidities greater than 
60 MV the limiting fluxes will be smaller. These re- 
sults are shown in Figure 3. Where only integral solar 
neutron intensities were measured, e.g., Hess and Kaifer 
(1967), the results have been converted to a differential 
neutron spectrum using the spectral shape presented by 
Lingenfelter and Ramaty (1967) and the known energy 
dependent response functions. At lower energies the 
OGO-6 result has the lowest experimental upper limit 
for the continuous emission of solar neutrons in the 
1-20 MeV neutron energy range. 

4. DISCUSSION OF RESULTS 

The OGO-6 solar neutron results place a lower upper 
limit in the 1-20 MeV energy interval and should be 
interpreted along with the recent higher energy results 
of Cortellessa et al . (1971), the limits set by Hei dbreder 


FIG. 3 


et al . (1970) and Eyles et al. (1971), and the earlier 
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measurements by Webber and Ormes (1967), all of which 

are shown in Figure 3. In the energy range 20-200 MeV 

the new limit set by Cortellessa et al . (1971) is within 

a factor of two of the upper limits previously set, 

which is very consistent agreement. In fact, no lower 

limits have been assigned in the 100-300 MeV energy 

interval than the early ones deduced by Webber and Ormes 

(1967). Although the calculations of Alsmiller and 

Boughner (1968) suggest that efficiency used by Webber 

and Ormes for the production of secondary protons by 

2 

solar neutrons in the first few grams/cm of the atmosphere 

might have been 50% too large, the value deduced for the 

solar neutron flux is still close to these new lower 

limits because the count rate limits set were very 

conservative. In the 2-20 MeV the OGO-6 solar neutron 

results set an order of magnitude lower upper limit than 

the measurements of Forrest an d Chupp (1969) and 

Coretllessa et al . (1971) in the overlapping energy range. 

The lowest measurements in Figure 3 are about 

comparable to the time-average intensity of the solar 

neutron flux from the many flares occurring during the 

last solar cycle 1954-1965 as calculated by Lingenfelter 

et al . (1965b). At 40 MeV this estimated flux was 

-5 2 

^>2 x 10 neutrons/cm sec MeV. Since the solar activity 
in the present cycle (1965 to ^1974) is considerably less 
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than the preceding one, the time-average intensity should 

be considerably less. This time-average intensity is 

still more than an order of magnitude greater than the 

theoretical upper limit obtained by Roelof (1966) based 

upon the IMP-1 protons, assuming all the protons resulted 

from neutron decays and the protons suffered isotropic 

diffusion. This suggests that the current measurements 

are predominately just background measurements. A 

comparison of the solar neutron flux to the atmospheric 

neutron flux supports this conclusion. At the geomagnetic 

equator for neutron energies E >50 MeV both the solar 

n 

and atmospheric neutron spectra have a similar energy 

dependence and magnitude if Lingenfelter and Flamm (1965) 

time-averaged values for solar neutrons are used. For 

5 < E <50 MeV the atmospheric neutron flux is two 
n 

orders of magnitude greater. Indeed, recent atmospheric 

neutron measurements by Simnett (1971) indicate a 

spectral shape which is relatively flat from 10 MeV to 

_ 2 

50 MeV than decreasing as E for E >50 MeV. This would 

n 

make it difficult to separate the two neutron sources by 
studies of the spectral shape at lower energies. 

More meaningful upper limits on the solar neutron 
flux can only be determined by neutron detectors with 
much larger geometrical factors, better efficiency and 
directionality to reduce the background of atmospheric 
neutrons. Large detecting systems which do not separate 
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neutron produced events from the background of gamma-ray 
interactions and other sources cannot provide more 
significant limits to the solar neutron flux unless 
detailed statistical analyses of the background counting 
rates are made to eliminate any possible fluctuations, 
as does in the analysis of Forrest and Chupp (1969) . 

Of course, these results do not exclude the possi- 
bility that there are neutrons from the sun during solar 
particle events. Unfortunately, no detecting systems 
with a high efficiency for neutrons have been aloft 
during any large solar flare events. In the 0G0-6 
measurements we looked for solar neutrons during the 
solar flare events between June and December, 1969. No 
positive evidence for any neutrons after several solar 
flares was found. No limits were placed upon the flux 
because the statistical uncertainty was too large. 

We conclude that the quiet-time solar neutron flux 

-3 2 

from 1-20 MeV cannot be greater than 1.8 x 10 n/cm 
sec at the 95% confidence level. This estimate is 
based upon a solar charged particle spectrum with a 
characteristic rigidity P = 


60 MV. 
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Figure 1 


Figure 2 


Figure 3 


FIGURE CAPTIONS 

- Neutron sensor arrangement and orientation 
on the 0G0-6 spacecraft. 

3 

He neutron counter efficiency for an isotropic 
neutron flux of energy. The details on these 
measurements are presented by Jenkins et al . 
(1970) . 

- Upper limits to the differential solar neutron 
flux at earth. Not all measured values are 
shown. For E >50 MeV the upper and lower 
limits are indicated with a band in which many 
results not cited lie. The smooth curve 
marked Lingenfelter and Flamm (1965) is the 
calculated time-average solar neutron flux at 


earth for solar cycle 19. 
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